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THE CMC FLUID MECHANICS COMPUTER PROGRAM SYSTEM IS BEING 
DEVELOPED TO TRANSMIT THE THEORETICAL EVOLUTION OF FINITE 
ELEMENT NUMERICAL SOLUTION METHODOLOGY# APPLIED TO NONLINEAR 
pi£LD PROBLEMS INTO A VERSATILE COMPUTER CODE FOR COMPREHENSIVE 
FLOWFIELD ANALYSIS. THIS REPORT IS VOLUME III OF A THREE 

VOLUME SET AND PRESENTS A DETAILED VIEW OF THE CODE FROM THE 
STANDPOINT OF A COMPUTER PROGRAMMER'S USE. 

A SYSTEM MACRO FLOW CHART IS PRESENTED ALONG WITH DETAILED 
FLOW CHARTS OF SEVERAL ROUTINES WHICH ARE DEEMED NECESSARY TO 
INTERACT WITH A THEORETICIAN / USER TO MODIFY THE OPERATION OF 

THE PROGRAM. 

DESCRIPTIONS OF ALL SUBROUTINES ARE GIVEN AND# WHERE PERTINENT* 
DETAILS OF USAGE# PRIMARILY FOR INPUT AND OUTPUT ROUTINES. 

DESCRIPTION OF INTEGER AND REAL SCALARS ARE GIVEN ALONG WITH 
A CROSS REFERENCE LIST DENOTING SUBROUTINE USAGE FOR THESE SCALARS. 

DESCRIPTION OF ENTRY POINTS IN DYNAMIC STORAGE VECTOR 'll'* 

THE LENGTHS OF EACH VECTOR ACCOMPANY THE SCALAR DEFINITIONS. 


A LISTING OF THE ROUTINES PECULIAR TO THE STANDARD TEST CASE 
IS INCLUDED ALONG WITH A LISTING OF THE INF'UT DECK AND PRINTOUT 
(WHICH INCLUDES DEBUG INFORMATION) FOR THIS CASE. 

THE THEORETICAL DESCRIPTION FOR 3DPNS IS GIVEN IN VOLUME I. 


THE USER ORIENTED DESCRIPTION IS GIVEN IN VOLUME II. 
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INTRODUCTION 


ALL VARIANTS OF THE CMC SYSTEM ARE BUILT UPON THE MACRO- 
ATAUCTURE ILLUSTRATED IN FIG* i. THE MAIN EXECUTIVE ROUTINE 
AU OCAiei CORE USING A VARIABLE DIMENSIONING SCHEME. BASEDURON 
?HE TOTAL DEGREES OF FREEDOM OF THE PROBLEM. THE SIZE OF THE 
i AfftcsT PROBLEM THAT CAN BE SOLVED IS THUS LIMITED BY THE CORE 
SIZE OF THE COMPUTER IN USE* THE PRECISE MIX BETWEEN NUMBER OF 
DEPENDENT VARIABLES (AND PARAMETERS) » AND FINENESS OF THE DIS- 
CRETIZATION IS USER-SPFCIFIABLE AND WIDELY VARIABLE* 

THE INF'UT MODULE SERVES ITS STANDARD FUNCTION FOR ALL 
DEPENDENT VARIABLE » PARAMETER AND GEOMETRIC COORDINATE ARRAYS* 

THE DISCRETIZATION MODULE FORMS THE FINITE ELEMENT 
DISCRETIZATION OF THE SOLUTION DOMAIN. 

THE INITIALIZATION MODULE COMPUTES THE REMAINING INITIAL 
PARAMETRIC DATA REQUIRED TO START THE SOLUTION. 

THE INTEGRATION MODULE CONST I TUTES THE PR I MARY EXECUTION 
SEQUENCE OF PROBLEM SOLUTION* IT IS BASED UFON AN IMFLICIT 
INTEGRATION ALGORITHM FOR THE COLUMN VECTOR OF UNKNOWNS FOR 

Solution which the discretized description is initial valued. 

CALLS TO AUXILIARY ROUTINES FOR PARAMETER EVALUATION# E.G.# 
VICQSITY# STRESSES# PHI# PP» ETC.# AS SPECIFIED FUNCTIONS OF 
DEPENDENT AND/OR INDEPENDENT VARIABLES ARE GOVERNED BY THE 
INTEGRATION MODULE* 

THE USER HAS CONSIDERABLE LATITUDE TO AD AFT 3DPNS ^OHOC 
TO THE SF'ECIFICS OF THEIR PARTICULAR PROBLEM AT THIS POINT BY 
DIRECTLY INSERTING USER WRITTEN SUBROUTINES INTO THE FLOW 
STREAM TO COMPUTE SPECIAL FORMS OF THESE PARAMETERS. 

THE OUTPUT MODULE IS SIMILARLY ADDRESSED FROM THE 
INTEGRATION SEQUENCE AND SERVES ITS STANDARD FUNCTION 
VIA A HIGHLY AUTOMATED ARRAY DISPLAY ALGORITHM. 
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PROGRAM DESCRIPTION 


THE 3DPNS VARIANT OF CMC IS WRITTEN TO BE READILY 
EXECUTED ON ANY LARGE SCALE COMPUTER SUCH AS THE IBM 370* 3031 f 
CYBER- 175* CYBER-203 OR UNIVAC 1108. IT IS PROGRAMMED IN 
FORTRAN IV EXCLUSIVELY * EXCEPT FOR CERTAIN MACHINE “DEPENDENT 
ROUTINES* E.G.* DATE* ADDRESS » ERROR HANDLING* ETC. 

THE PROGRAM USES ABOUT 100K WORDS OF MEMORY FOR A SOLUTION 
REGION INVOLVING ABOUT 200 NODES AND ABOUT 12 DEGREES OF 
FREEDOM PER NODE. OUTPUT CAN BE STORED ON DISC* TAPE OR PUNCHED 
CARDS FOR FUTURE RESTARTS. 

CONSIDERABLE EFFORT HAS BEEN DEVOTED TO CONSTRUCTION OF 
CMC/3DPNS IN A UNIFORM AND CONSISTENT MANNER. 

A COMPREHENSIVE SUBROUTINE SUBSTRUCTURE IS UTILIZED 
TO ALLOW FUNCTIONAL PROCESSING IN A CONSISTENT AND READILY 
FOLLOWED SEQUENCE. INPUT FORMATS ARE CONSISTENT AND 
FLEXIBLE* AND PROGRAM INPUT CONTROLS ARE STANDARDIZED. 
CONSISTENT NOTATION OF VARIABLE NAMES IS EMPLOYED IN ALL 
SUBROUTINES AND ONLY THREE MAJOR ARRAYS ARE MAINTAINED. 

THE COMMON BLOCK / VARBLE / CONTAINS THE ' I ARRAY ' AND 
THE 'RARRAY' VECTORS OF SCALARS. THE I ARRAY CONTAINS A LIST 
OF ALL THE INTEGER SCALARS USED IN COMMON THROUGHOUT THE 
PROGRAM. THE RARRAY CONTAINS A LIST OF ALL THE REAL SCALARS 
USED IN COMMON THROUGHOUT THE PROGRAM. 

THE COMMON BLOCK / ARRAYS / CONTAINS THE 'IZ' ARRAY. 

THIS ARRAY CONTAINS ALL THE VARIABLE DIMENSIONED VECTORS* 
MATRICES AND HOLLERITH DATA WHICH ARE USED BY THE PROGRAM. 

THE SOURCE LISTING OF CMC/3DPNS IS NOT INCLUDED IN THIS 
DOCUMENT DUE TO ITS EXCESSIVE LENGTH < APPROX. 17000 CARDS >. 
THE SUBROUTINES FOR THE 3DPNS VARIANT OF THE COMOC PROGRAM ARE 
DESCRIBED IN THIS MANUAL* 
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THE FOLLOWING IS A LIST OF COHHON BLOCKS WHICH ARE USED 
TO STORE DATA WHICH HAY BE DESTROYED BY OVERLAY PROGRAMS t 
OR WHICH IS NEEDED ACROSS CERTAIN SUBROUTINES* 

THE LENGTH OF THE COHHON BLOCKS ARE SHOWN IN PARENTHESIS* 


/ARRAYS/ (SET IN HAIN) 

/CONPPR/ <10> BDINPTt BLOCK DATA* CONTES t DERVDXt XHPLCT t IHPSLV t 
PPRES» PRSGRDr RNLDST t TCHECKt 

/DERIV / (100) BCONDT t DERVBL t DFCFBLt DRVBUGt IHPLCTt IHPSLV t 
INDEX t PPRESt STCODEt 
/DISCRT/ (10) BLOCK DATA* DSCRTZt 
/IHPDAT/ (20> BDINPTt HNK2t IHPSLV t 
/JADRES/ (30) 

/NEWGOH/ (10) BLOCK DATA* NWGEOH » 

/NPINTQ/ (30) 

/OUTDUH/ (75) BDINPTt BLOCK DAT At REOUTPt 
/PLOTFE/ (46) BDINPTt BLOCK DATAt FEPLOTt 
/PPRSET/ (41) BDINPTt LlNK2t GETPPRt 
/PRGDUH/ (50) PRSGRDt 

/RDREND/ (1) BDINPTt ADDDELt DELELHt REDRELt 

/STDHAT/ (30) BDINPTt ACONDT t DERVBL t FEDIHNt PPRESt STRFt 

/VARBLE/ (1000) 


C ■ • t* 

C:1 «Ll ^ 

DESCRIPTION OF CONTROL WORD INPUT (BDINPT) 


A CONTROL CARD WITH THE FOLLOWING PARAMETERS IS READ IN - 


PARAMETER 

FORMAT 

CARD COLS 

. DESCRIPTION 

VI 

A8 

1-8 

CONTROL VARIABLE* 

NMUL 

FREE 

AFTER 8 

Nl* N2» . . ♦ * NN 


NMUL IS A VECTOR 

'N' LONG* IT CAN CONTAIN El THE! 


| INTEGER OR REAL VALUES. 


FOR CLARITY OF READING THIS SECTION* THE FOLLOWING 
EQUIVALENCES WILL BE ASSUMED IN SUBROUTINE BDINPT ♦ 


FI* 

N1 

8 

F2» 

N2 

ss 

F3» 

N3 

s 

F4* 

N4 

s 

F5» 

NS 

s 

FA* 

NA 

e 

F7* 

N7 

8 

F8» 

N8 

s 

F9* 

N9 

s 


FMUL<1) » 
FMUL<2) » 
FMUL < 3 ) * 
FMUL<4 ) » 
FMUL(S) » 
FMUL < 6) * 
FMUL (7) » 
FMUL (8) * 
FMUL(9) * 


NMUL(t) 
NMUL<2> 
NMUL ( 3 ) 
NMULC 4 > 
NMUL ( 5 > 
NMUL (6) 
NMUL < 7 ) 
NMUL (8) 
NMUL (9) 


IF KDUMP = 1 IN NAMEOl* THEN THE ENTIRE INPUT CARD IS PRINTED 
IMMEDIATELY AFTER BEING READ FOLLOWED BY THE DATA THAT IS BEING 
STORED ALONG WITH THE DATA'S ENTRY POSITION IN THE 12 ARRAY 
(SEE IAINTG AND I AREAL) 


THIS ROUTINE LOOKS FIRST FOR A MATCH OF VI WITH CERTAIN KEY WORDS 
WHICH WILL EITHER CAUSE A SUBROUTINE TO BE CALLED OR PROGRAM 
FLOW TO OCCUR* 
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DESCRIPTION OF CONTROL WORD INPUT (BDINPT) 


THE KEY WORDS THAT ARC SCANNED ARE - 


1 . 

2 . 

3* 

4. 


6 ♦ 
7. 
8 * 
9. 


(BLANK) 
COMMENT 
COMOC 
COMTITLE 
DESCRIPT N1 


RETURN TO SCAN ANOTHER CARD. 

PRINT THIS CARD AS A COMMENT. 

PRINT THE COMOC TITLE PAGE TWO TIMES. 

READ A TITLE CARD WHICH WILL APPEAR ON COMOC. 
CALL DSCRPT AND PROCESS ACCORDING TO N1 . 


END 

EXIT 

FEDIMN 

FENAME 


N1 

BLANK 

203 

204 


READ AND WRITE INFORMATION CARDS. 

READ TITLES FOR DEP. VAR. OUTPUT HEADING. 
READ DESCRIPTIVE TITLE FOR HEADING 
AT BEGINNING OF OUTPUT HEADER. 

332 - READ TITLES FOR PARAMETERS PRINTED IN 
THE OUTPUT HEADER. 

- RETURN CONTROL TO MAIN PROGRAM. 


- CALL EXIT. 

- CALL DIMENSIONALIZATION ROUTINE FEDIMN. 

- CALL FENAME TO SET DEFAULT SCALARS AND THEN 
CALL NMELST TO READ IN NAME01 AND NAME02 


10 * 

11 * 

12 . 

13. 


14. 

15. 
16* 

17. 

18. 

19. 

20 . 

21 . 

22 . 

23* 

24. 

25. 
26* 


ICOND 

I ARRAY Nl» 

INPUT N1 
KBNO N1 

KBNO N1 t 1 


NAMELISTS. _ _ A| A __ 

- CALL ICOND TO PRINT INTEGER AND REAL SCALARS. 

N2» N3» N4» ETC. 

- SET IARRAY(Nl) * N2» I ARRAY <N3> = N4> ETC. 

- SET INPUT UNIT TO HI. 

- ENTER FIXEP NODES FOR DEP. VARIABLE N1 
BY CALLING SUBROUTINE GETBND. 

- ENTER FIXED NODES FOR DEP. VARIABLE Nl 

BY CALLING SUBROUTINE GETBCD WHEN N2 .EQ. 1. 


LINK1 N1 
LINK2 N1 
LINK3 N1 

LINK4 N1 
LINKS N1 
NAMELIST 

PLUS N1 N2 N3 

QKNINT 

RARRAY N1 N2 . 

VYYEND N1 
ZZZZZZ 
2DBL 
2DPNS 


- CALL LINKKN1) 

- CALL LINK2(N1 ) 

- CALL LINK3(N1 t U, Y> WHERE U» Y HAVE MEANING. 

WHEN N1 » 1. k 

- CALL LINK4(N1 » K) WHERE K IS A DUMMY INTEGER. 

- CALL LINKS (Nl) 

- CALL ENTRY ' NMELST ' IN SUBROUTINE FENAME 
TO READ IN NAME01 AND NAME02 NAMELISTS. 

. . . ETC. 

R(N1) ■ R(N2) + R(N3) + . ♦ . + R(NN> 

START THE INTEGRATION PROCESS. 

. . Ft . ♦ ♦ rN ♦ ♦ . NN 

- R(N1> * R(N2> * .FI * * FN * . # R(NN) 

- DENOTES END OF INPUT FOR DEP. VAR. Nl. 

- CALL SUBROUTINE ZZZZZZ. 

- SET NM * 2. 

- SET NM » 2. 
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DESCRIPTION Of CONTROL WORD INPUT 


\ S S \ VECTOR OPERA I IONS \ \ N 

THE FOLLOWING INDICES WILL PREVAIL WHEN SUBSCRIPT NOTATION IS USED, 


IZ1 

8 

IZ< N2 ) 

+ 

N6 

IZ2 

S 

IZ(N3) 

+ 

N7 

IZ3 

8 

IZ<N4) 

+ 

N8 

IZ4 

8 

IZ(N3) 

+ 

N4 

IZS 

8 

IZ(N2> 

+ 

N3 

IZ6 

8 

IZ(N4> 

+ 

NS 

IZ7 

8 

IZ(N6) 

+ 

N7 

IZ8 

8 

iz<n:o 

+ 

N4 

IZ9 

8 

IZ(N3) 

+ 

NS 

R(N) 

8 

RARRAY(N) 


I (N) 

8 

IARRAY(N) 



£25??2 D NAME IS LOCATED IN COLUMNS 1 - 8 STARTING IN COL. 1 
ENTRIES START AFTER COLUMN 8. COL * 14 


COMMAND [ ENTRIES 


NAME 


DESCRIPTION 


1, ABSVAL N1 N2 N3- SET RZ(I2<N2>) = ABS<RZ< 1Z<N3> > . I « 1> N1 

2. INDEXSET N1 N2 N3 N4 NS N6 N7 N1 

** DO FOR I * 1, N1 
K = IZ(IZ5) - 1 
RZ( IZ6+K) * RZ ( IZ7+I > 

ALLOWS NEW VALUES TO READ INTO A SEQUENCE 
VECTOR ATIONS IN ™ E VINTeR ' IZINTG AND IAINTG 

CALL MATSUM ( RZ(IZ(N2)>» RZ(IZ(N3))» R<N4>» 
RZ< IZ<N5) > » N1 ) 

CALL PDUMP ( IZ(IZ(N2))r IZ(IZ(N3))r N1 ) 


3* INTEGER 


4. MATSUM 


5. PDUMP N1 N2 N3 - 

6. READ N1 N2 N3 N4 


7. REAL 


IF ( N2 .EQ. 0 > N2 * NNODE 
IF ( N2 .GT. 0 ) N2 ■ IARRAY(N2> 

IF ( N4 *EQ* 999 ) IZ4 * N3 
READ < Nl ) RZ( IZ4) r I * 1» N2 
- ALLOWS NEW VALUES TO READ IN A SEQUENCE 

2!Li 0 £ ATI0N8 IN THE IZREAL AND IZREAL 

VECTORS. 

8* RECIP Nl N2 0 0 0 N6 

CALL RECIP ( Nl. RZ<IZl> > 

10* REBTOR Nl N io N2 ” £?!“!“ 2?™? <lf IZ<N2> ' «<ITY) ) WITH NP - Nl 
11. SETVAL Nl N2 N3 N4 NS ' 

CALL SETVAL < RZ(IZl>r RZ(IZ2)> R(N4>» R(NS>» Nl ) 
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dlocriptiom or comtroi word input 


i?., s;qrt iu N2 ... N.J ... NO 

R(Ni) R<N2>mXPNI 
R(N3> R<N4)*#I:XPNT 

. . . me) dy paifus. <i. : vem mo. oi y h i r : i o > , 

FOR AN ODD NUMBER 01 ENTRIKO* THEN 

SORT Mi Mi* M3 - F<2U20> « R2< 129) **EXPN I . I > Mi 

WHERE EXPNT “ 0.1) UNLESS 'ROWER FI' CARO F*RE CEDES 

I HE SORT CARO* (HEM liXPNT II. 

13. VECMUL Ml Mi* M3 N4 

RZCm> « R2<1Z2> # R2<I23)**I TOWER 

WHERE F POWER « 1.0 OR F9* IF FV .ME. 0.0. 

14. WRITE Ml M2 M3 N4 

IF < M2 .EU. 0 ) M2 « NNODE 

IP < M2 .OF. 0 ) M2 « IARRAY(M2> 

IF ( M3 ,0T, 0 < AMU . Ml «E0. A ) 

CALL 001 VE C < M2 1 17(124 >, PIE ) 

IF ( M3 .LI. 0 .AND, Ml ,E0, A > 

CALL 00 1 MOD < N2r 12 < 124)* P I I- > 

IF < Ml ,E«. 7 > 

WRITE (7,911*40) ( (IU(I24 H -1 > » I« 1 rN2) 

IF ( Mi .ME, A .AND, Mi .ME, 7 ) 

WRITE (Mi) ( <R2( I24-H-1 ) r I< ! 1 » M2) 


DESCRIPT ION OK INTEGER TYPE INPUT < RDJNPI > 

THE FOLLOWING COMMANDS ARE UTILIZED TO LOCATE ~STOR AGE ~ADDRE ssp *i ~ 
FOR INTEGER ARRAY DATA SPECIFICATION. THE CARDS ARF DIREHLY E 
FOLLOWED BY INTEGER DATA CARDS IN FREE FORMAT* THE READ SCAN IS 
TERMINATED BY A ' T ' OR BLANK CARD. 8tAM lb 

* * * NOTE * * * 

■*' ?he°mm M Ihd assr uhich cum,mns 

thk i? ftRRftv ftT WH,CH 10 •*“» 

IARRAY<XAlNTG(IO) «• NUMBER OF ENTRIES STORED STARTING AT 

Jr [I? * EQ * ENTER LITERAL DATA, (SEE SUBROUTINE BETBND) 

IF N1 «NE. 1 t ENTER NUMERICAL DATA DIRECTLY. 

EXAMPLE? 

I OS AYE 

1248 2248 3248 1247 9248 * 5248 6248 7248 T 

IOMULT 

3*2 21 4*2 » 8*1 T 


K 

COMMAND 

I Z INTO 

IAINT 

mm _ 

NAME 

<K> 

(K> 

7 

IPINT 

5 

31 

12 

LINKCALL 

121 

125 

14 

IOMULT 

123 

67 

15 

IOSAVE 

124 

60 

16 

CNTPTS 

127 

47 


DESCRIPTION 


17 

CNTNDS 

128 

128 

18 

XBORD 

38 

131 

19 

XPSXBD 

14 

67 

21 

IONUMB 

131 

142 

22 

MPARA 

135 

67 

28 

ELEMENTS 

26 

67 


SOLUTION SEQUENCE VECTOR. 

LINK NOS. TO BE CALLED AT END OF IMPLCT 
OUTPUT VARIABLE MULTIPLIER FROM R ARRAY, 
VARIABLE LIST TO BE DISPLAYED AT OUTPUT 
CONTOUR NODES T. BE USED IN 

CONTES r DFCFB * TRBTHKr WLFLXSr ETC. 

* * * NOTE * * * 

WHEN CNTPTS AND CNTNDS ARE USED* IARRAY<47) 
SHOULD BE SET TO AT LEAST THE NO. OF ENTRIES 
IN CNTPTS TO ALLOW FOR STORAGE ALLOCATION. 
NO. OF NODES IN EACH CONTOUR LINE. 
COUNTER-CLOCKWISE LIST OF BOUNDARY. 
DIAGONAL FLAG REVERSAL VECTOR. 

LIST OF ENTRIES IN RARRAY TO BE 
DISPLAYED AT START OF EACH OUTPUT. 
LIST OF MULTIPLIERS IN RARRAY USED 
TO MULTIPLY IONUKB ENTRIES. 

READ IN ELEMENT NODE CONNECTIONS. 


* * * NOTE - 

IAINTG(K) OR IAREAL(K) 


67 IMPLIES PRESET LENGTH IS NOT CHANCED. 


t)F POOR QUAUTY 
JEgCRIPTl OW or REfl l- TTFj: iNr ui (BBINPT) 


FOR REAL^RRAY^ATA^^^ NAME ^Aer? 1 ^ AW,RESSE S 

FOLLOUEO BY REAL BATA CARDS IN FREFFOr 22? ?SS® e ?^ " IKEC11 ' 
TERMINATED BY A *T' OR BLANK CARI.7 T- HE REftr ' SC *N IS 


* * * NOTE * * * 

THE 0 MMmInd NAMES I FOR I REAL E TYFE T ?NF'UT? EftL ’ UHICH C0NTAINS 
,Z, 8 ?Sr?N0 K REAl‘ ENTRIES** ,M ™ E 12 AR ^ Y AT “ HICH T ° »«• 
IAR IZ(i1rEAL(K))! ' NUM8ER 0F ENTRIES ST0RE * STARTING AT 

™e S en?ereS°da?a? DS CAN contain " 0R0UR ° F "»-™ ‘RS T.,< 

THE i? l I/? l rml e !1J? IE8 can “ 0F Th ***•» 

3 2 : -n n : titles rarray(N) :i k: er - 

EXAMPLE! 

DEPVAR 2 3 iOO.O -27 

30*22 ♦ S 2 6.7 35*3 10*50*0 T 

MULTIPLIES EACH VECTOR ENTRY OF HEP. VAR. **> by 
100.0 * RARRAY ( 3 > / RARRAm?) 


KjCOMMAND flZREAL ^ IAREAlI DESCRIPTION. 

1^ VU2P0S I 

2 l VU2VAL I 

3* UU3P0S | 

I 

4' UU3UAL { 

9 i UX1C0R J 
10 « VX2C0R j 


(K) 

| CK) 

63 

f 177 

6A 

f 67 

65 

| 178 

66 

1 47 

89 

' 16 

90 1 

1 16 


r X F0R VARIABLE GRID CHANGE 

f erI^ T 2 ANSVERSE direction. 

SCALE FACTOR FOR GRID CHANGF 

I v I? ™5" SVERSE DIRECT*™. 


... ' wn vrmAHUUI 

IN NORMAL DIRECTION! 
:ale factor for mark 

IN NORMAL DIRECTION, 
DEFAULT « 1.0 / ALC 
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DESCRIPTION OF REAL TYPE INPUT 


K I COMMAND flZREAL I 1 AREAl! DESCRIPTION* 


mm m 

NAME j 

<K> | 

<K> 

16 | VPRESS 1 

91 1 

67 

18 [ 

VYY 1 

72 1 

1 67 


VSUTHLD , 

133 

67 

( 

( 

f 

25 

j 

VX3ST j 

139 ! 

I 161 

1 

26 

VPVSX | 

140 , 

1 

1 47 

29 1 

RARRAY 

0 

I " 

33 | 

| 

DEPVAR 

• 

NN 

i 47 


I PRESSURE VALUES AT NODE POINTS. 
t DEFAULT * PINF 

[DEPENDENT VAR. DIST. AT NODE POINTS. 
STLDVR t STLDTR* STLDCRt STLDEXt STLCON 
| ENTRIES FOR SUTHERLANDS LAW. 

DEF. . 1163E-4» 494.0* 204.0» 1.5* 0.0 
! DOWNSTREAM STATIONS AT tlHICH PRESSURE 
IS DEFINED. 

[DOWNSTREAM PRESSURES AT VX3ST » 

DEFAULT « XMUINF 
JRARRAY(NX) * AMULT t WHERE 
'aMULT * COMBINA*. ION OF REMAINING ENTRIES. 
.DEPENDENT VAR. DIST. AT NODE POINTS. 

Inn is determined by def*. var. ni. 

iNN = IYY + (NP-1 > # NYY * NODE 
1 WHERE NP IS THE LOCATION OF Nl IN 
I THE IPINT VECTOR. 

THE AMULT MULTIPLIER FOR THIS CASE CONSISTS 
OF THE PRODUCT OF THE REMAINING CONTROL 
NOS* AFTER Nl. 


ORIGINAL Pi’.va 13 
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SUBROUTINE DESCRIPTIONS. 


THE FOLLOWING PAGES CONTAIN A BRIEF DESCRIPTION OF 
THE SUBROUTINES IN THE COMOC COMPUTER PROGRAM. 


NAMES IN PARENTHESES INDICATE CALLING ROUTINES. 


IF NO NAME IS ENTERED THEN SEVERAL ROUTINES PLACE CALL 


-f’.VCIS 

OF Huoft QUALITY 


CONTROL ROUTINES 


-MAIN 

THIS IS THE MAIN CONTROL PROGRAM WHICH INITIALIZES THE RARRAY, 
I ARRAY AND THE IZ ARRAYS TO ZERO. 

TO CHANGE THE VARIABLE SOTRAGE CAPACITY OF THE IZ ARRAY* RESET 

I25frF X ?5?? ION ° F 12 AND ' ACCORDINGLY, THE VALUE OF 'I7SIZF'. 
AFTER INITIALIZATION THE CONTROL ROUTINE BDINPT IS CALLED. 


PPINPT (MAIN) 

THIS IS THE CONTROL ROUTINE WHICH INITIALIZES VECTORS AND 
CONTROLS THE FLOW OF THE PROGRAM ACCORDING TO USER INPUT. 
( SEE TABLE 1 FOR INPUT DESCRIPTION. ) 


j- INK1 SUB. LINK1 (KK) 

PLACE CALLS TO THE FOLLOWING ROUTINES 
ACCORDING TO THE VALUE OF KKJ 

2. NODELM 

3. GEOMFL 
A. GETPPR 
5. PRSGRD 
7. BCONDT 


•LINK? SUB. LINKS (KK) 

PLACE CALLS TO THE FOLLOWING ROUTINES 
ACCORDING TO THE VALUF OF KKJ 

1. DFCFNS 

2. DFCFBL 

3. WLFLXS 

4. CONTES 

5. REOUTP 

A. TRANSFER TO FEOUTP PORTION OF REOUTP. 
7. STRF 
9* DRHOBL 
10. TBLINP 
13. XYCRDM 
14* DSCRTZ 
15. TRBT! ; 

23. CPSTUP 


n 


▼ 



PAGE IS 

OF POOR QUALITY 


CONTROL ROUTINES 


LINK3 SUB . LINN3 ( KK» U> Y ) 

PLACE CALLS TO THE FOLLOWING ROUTINES 
ACCORDING TO THF VALUE OF KKJ 

1. NBNDRY 

2. RITE 

4. DIMEN 

5. SUMKEY 


LINK4 SUB. LINK4 < KK» NN ) 

PLACE CALLS TO THE FOLLOWING ROUTINES 
ACCORDING TO THE VALUE OF KKJ 
2. 1MPLCT 
9. POTENT 


LINKS SUB. LINKS ( KK ) 

PLACE CALLS TO THE FOLLOWING ROUTINES 
ACCORDING TO THE VALUE OF KK J 
1* NUGEOM 

2. CALL ROUTINES FROM LINKCALL LIST AT END OF IMPLCT. 

3. CPINIT 
5. SCHPRN 
A. SETDIF 


k 

( 




> 

♦ 

f 


I 


i 
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•V 


ORtQW.*\ p,«\c2 jg 

OF POOfJ guALITY 


GENERAL PURPOSE ROUTINES 


BOOLEAN ASSEMBLY OF ELEMENT MATRIX INTO GLOBAL MATRIX* 

AVR6 SUB. *VRG ( IND* ARR* NUMB > 

COMPUTE THE ARITHMETIC AVERAGE OF 'NUMB' ENTRIES 

IN THE ARRAY 'ARR' INDEXED BY THE ENTRIES OF 'IND'* 

-CSFINP (NODPCP) 

CUBIC SPLINE FIT FOR A MONOTONICALLY INCREASING FUNCTION. 

DUDY (ULFLXS) 

THREE POINT INTEGRATION FORMULA FOR COMPUTING DUDY. 

LOC (MPRD) 

COMPUTE VECTOR SUBSCRIPT FOR AN ELEMENT IN A MATRIX OF 
SPECIFIED STORAGE MODE. 

LOCATE SUB. LOCATE ( NA» NN* M* N ) 

FIND THE LOCATION OF 'M' IN THE ARRAY 'NA' AND STORE IT IN 'N'. 

LOOK SUB. LOOK < M* TIN* CIN* TFIN* COUT > 

FIND 'TFIN' IN 'TIN' VECTOR. PERFORM LINEAR INTERPOLATION 
OF 'TIN' * 'CIN' AND PLACE IN 'COUT'. 

LOOKAV < NUGEOM ) SUB. LOOKAV ( M» TIN* CIN* TFIN* COUT* OSNEL. > 

SIMILAR TO 'LOOK' EXCEPT THAT TRANSITION POINTS ARE RELAXED 
BY SCALAR 'RELAX' TO PROVIDE 'SMOOTH' TRANSITION OF GEOMETRY. 

MATSUM SUB. MATSUM ( A* B* COEF* C» H ) 

I ARRAY (347) = INCLFT 

I ARRAY (348) * INCRGT 

KRINC = 1 

NN - N * INCLFT 

DO 200 I « 1* NN* INCLFT 

A( I > > B( I ) 4 COEF * C(KRINC) 

KRINC - KRINC + INCRGT 
200 CONTINUE 

INCLFT ■ 0 
INCRGT - 0 




OF PU'.'li QUALIFY 


GENERAL PURPOSE ROUTINES 


MINMAX ( PSCRTZ » ORDER » SETSCL ) 

SUB. HINMAX (XN* XM* ILIM* NAR* INN. IMX > 
COMPUTE THE MINIMUM 'XN' AND MAXIMUM 'XM' ENTRIES IN AN ARRAY 
' NAR ' AT LOCATIONS 'IMN' AND 'IMX' IN THE ARRAY. 

UNMX (ELEM) SUB. MNMX ( NN* INA* MN* MX* LOW* LHI ) 

FROM AN INTEGER VECTOR INA CONTAINING NN ENTRIES* 

STORE THE FOLLOWING: 

LOW - POSITION IN 'INA" OF MINIMUM. 

LHI - POSITION IN 'INA' OF MAXIMUM. 

MN - MINIMUM VALUE IN 'INA'. 

MX - MAXIMUM VALUE IN 'INA' 

MPRD (GEOMFL) 

MULTIPLY TWO MATRICES AND STORE IN RESULTANT MATRIX. 

C = A * B 

■NBNDRY (LINK3( 1 ) > 

THE VALUE OF ' NBSET ' DETERMINES THE OPERATION OF THIS ROUTINE* 

NBSET = 1* STORE THE DEPENDENT VARIABLE INTO 'YY' VECTOR. 

NBSET = 0* STORE THE DEPENDENT VARIABLE FROM 'YY' VECTOR. 

NODFIX <BDINPT *DFCFBL) SUB. NODFIX ( IB* INFIX ) 

COMPUTE FIXED NODES AS A FUNCTION OF INPUT IN BDINPT . 

IN DFCFBL* COMPUTE FIXED NODES FOR TKE - DISS. AS FCT . OF YPLUS. 

ORDER < COLS * ROWS »XYSCAL) 

“ORDER 3 ARRAYS ACCORDING TO THEIR XI AND X2 COORDINATES. 

THE THIRD ARRAY WILL CONTAIN THEIR ARRAY LOCATIONS. 

RECIP SUB. RECIP ( NN* A ) 

FOR NON-ZERO ENTRIES IN 'A'* STORE 1.0/A BACK INTO 'A'. 

RESET SUB. RESET ( NN* A* V ) 

RESET 'NN' ENTRIES OF ARRAY 'A' TO THE VALUE 'V'. 

RESET I SUB. RESET I ( NN* NVEC* NVAL ) 

RESET 'NN' ENTRIES OF ARRAY 'NVEC' TO THE VALUE 'NVAL'. 

BESTOR SUB. RESTOR ( NN* A* YY ) 

SETS DEP. VARIABLE 'NN' IF SOME ENTRIES ARE CHANGED 
WITHOUT INTEGRATION OR ITERATION. 

NI * FCT (NN) 

YY (NI+I ) * A(I) FOR I « 1* NNODE 


.A 




GENERAL PURPOSE ROUTINES 


3ETRZS SUB SETRZS ( IOUT, IN» C, D» N > 

COMPUTE RZ( IOL . tl-l > = C * RZUN+I-1) + D, FOR I = 1, N 

SETVAL SUB, SETVAL < A, Bt Ct D» N ) 

I ARRAY (347) * INCLFT 
I ARRAY (348) = INCRGT 
KRINC * 1 
NN * N * INCLFT 
DO 200 I * It NNf INCLFT 
A(I) « C * B(KRINC) + D 
KRINC = KRINC + INCRGT 
200 CONTINUE 

INCLFT = 0 
INCRGT = 0 


■M6BHAX (FEDIMN t FEPLOT ) SUB. VARMAX ( 

FOR +NK» FIND MAXIMUM VALUE IN VECTOR 
FOR -NK» FIND MINIMUM VALUE IN VECTOR 


NK» VAR> XMAX ) 
'VAR' AND STORE IN 
'VAR' AND STORE IN 


'XMAX 

'XMAX 


yECMUL SUB. 

FOR Iaj, NPT» STORE! 


VECMUL ( NPT» 


C(I) * A(I) * B<I> 


C» 


A * B 


) 


INPUT ROUTINES 


OR!C*' R -’" !S 

OF POOR QUALITY 


Mim, <ELEM,GETBCD, GETBND) 

VALUE%f E ^KTYP% TRIES IN AN INTE6ER ARRAY DEPENDING on the 

KTYP * 1, DELETE 
KTYP * 2, ADD 


BNDSET (GETBCD, GETBND) 

DETERMINE NODES TO BE INSERTED INTO BOUNDARY ARRAY. 


DESCRP<KA> (BDINPT) 
FOR DIFFERENT VALUES 
‘READ/WRITE 
READ/UR I TE 


0,1 

3 

203 

204 


READ TITLE 
READ TITLE 


SUB. DESCRP(KA) 
ka> rr *nt THE FOLLOWING. 
PROBLEM DESCRIPTION. 

I ARRAY ENTRY AND DESCRIPTION 
RARRAY ENTRY AND DESCRIPTION 
FOR OUTPUT STATION HEADING. 
FOR 'COMOC' COVER PAGE. 


( 10A4 , 110, SA4 l 
( 10A4, 110, SA4! 


LET BCD (BDINPT) 

INPUT Jn E s ^ck*pe“ ssr except that ™ ere ,s ° nly 


SAME DEF. AS IN GETBND. 
DESCRIPTION 


KEYWORD 

FORMAT 

BLOCK COLS. 

(SAME) 

A8 

1-8 

KEYWORD 

FORMAT 

BLOCK COLS. 

K0DE1 

FREE 

AFTER 8 

K0DE2 

FREE 

AFTER K0DE1 

K0DE3 

FREE 

AFTER K0DE2 

Al 

FREE 

AFTER K0DE3 

MAI 

FREE 

AFTER Al 

A3 

FREE 

AFTER MAI 

MA3 

FREE 

AFTER A3 

+MA1 * 
-MAI » 


SAME AS K0DE1 IN GETBND. 

SAME AS K0DE2 IN GETBND, 

SAME AS K0DE3 IN GETBND. 

JfJJ-UE 0F A1 F0R THIS boundary. 

RARRAY MULT. FOR Al. 

VALUE FOR A3 BOUNDARY CONDITION* 
RARRAY MULT. FOR A3. 

+MA1 * MULTIPLY BY RARRAY (MAI) 

DIVIDE BY RARRAY (MAI ) 


£ETBCH (GET BCD) 

EXTRACT BOUNDARY CONDITION VECTORS FROM INPUT DATA. 


INPUT ROUTINES 


on:wn r.vs js 

OF POOR QUALITY 


6ETBND (BDINPT) 

ESTABLISHES THE BOUNDARY NODE VECTOR FOR EACH DEP . VAR. USING 
EITHER THE WORD 'ADD' OR SIMPLE GEOMETRY OF THE PROBLEM 
WITH THE FOLLOWING KEYWORDS AND CODES - 
EACH CARD IS DIVIDED INTO FOUR IDENTICAL 
BLOCKS OF 20 COLUMNS EACH. 

ALL BLOCKS ARE OF THE SAME FROMAT SO THAT A 
DESCRIPTION OF ONE BLOCK ONLY WILL BE GIVEN. 

THE BLOCKS START IN COL. 1* 21 * 41 AND 61. 



KEYWORD 

FORMAT 

[block COLS. 

DESCRIPTION 


TOP 

A8 

1 - 

8 

ACROSS TOP FROM LEFT TO RIGHT. 


-TOP 

A8 

1 - 

8 

ACROSS TOP FROM RIGHT TO LEFT. 


BOTTOM 

A8 

1 - 

8 

ACROSS BOTTOM FROM LEFT TO RIGHT. 


-BOTTOM 

A8 

1 - 

8 

ACROSS BOTTOM FROM RIGHT TO LEFT. 


RIGHT 

A8 

1 - 

8 

UP RIGHT HAND SIDE. 


-RIGHT 

A8 

1 - 

8 

DOWN RIGHT HAND SIDE. 


LEFT 

A8 

1 - 

8 

UP LEFT HAND SIDE. 


-LEFT 

A8 

1 - 

8 

DOWN LEFT HAND SIDE. 


(BLANK) 

A8 

1 - 

8 

IGNORE BLOCK. 


ADD 

A8 

1 - 

8 

CALL ADDDEL TO INSERT ENTRIES. 
IGNORE BLOCK COLS. 9-20. 


DELETE 

A8 

1 - 

8 

CALL ADDDEL TO DELETE ENTRIES. 
IGNORE BLOCK COLS. 9-20. 


DONE 

A8 

1 - 

8 

LEAVE ROUTINE. 


FOR THE FOLLOWING KEYWORDS* THE THREE CODES (WE'LL CALL THEM 
K0DE1 * K0DE2 AND K0DE3 FOR CONVENIENCE) WILL DETERMINE WHICH 
NODES WILL BE SELECTED. 


KEYWORD 

FORMAT 

BLOCK COLS. 

DESCRIPTION 

K0DE1 

14 

9-12 

ROW OR COLUMN DISPLACEMENT FROM 




EDGE BEING DESCRIBED (DEF* = 0). 

K0DE2 

14 

13 - 16 

POS. IN LINE TO START (DEF. = FIRST) 

K0DE3 

14 

17 - 20 

POS. IN LINE TO END (DEF. ■ LAST) 


( BDINPT *DSCRTZ*READVJ SUB. READER ( LI* FI* KNTR > 

WHEN DATA CARD IS READ. STORE 

1. INTEGERS IN 'Ll' VECTOR 

2. REALS IN 'FI' VECTOR 


4 n 


ORIGIN - Pr.'-':E IS 

OF POOR QUALITY 


INITIALIZATION ROUTINES 


BRDSHW ( BDINPT ) 

DISTRIBUTE U1 PROFILE OVER DOMAIN ACCOR NO TO BRADSHAW DATA. 

- CASE NO. 1400 FROM 1967 STANFORD PROCEEDINGS . 

CPINIT < DIMEN) 

COMPUTE CPINF AT TSINF. 

.DIMEN (LINK3<4) ) 

COMPUTE NON-DIMENSIONAL I ZING FACTORS USED IN PROGRAM. 

FEDIMN (BDINPT) 

SET UP DIMENSIONS OF VARIABLE LENGTH ARRAYS USED IN THE SYSTEM. 
IF ' KDUMP ' * 1» PRINT LOCATION OF ENTRY POINTS IN 'IZ' ARRAY. 

FENAHE (BDINPT) 

THIS ROUTINE CONTAINS A LIST OF ALL EQUIVALENCED VARIABLES 
IN THE I ARRAY AND RARRAY VECTORS. 

MOST DEFAULT VALUES ARE ALSO SET IN THIS ROUTINE. 

CALL NMELST TO READ IN NAME01 AND NAME02 NAMELISTS. 

f 

FINDBE (BDINPT) 

DETERMINE A SERIES OF BOUNDARY ELEMENTS AS A FUNCTION OF 
INPUTING BOUNDARY NODES IN COUNTER-CLOCKWISE ORDER. 

ON FIRST PASSt IF IBORD IS READ* FIND BORDER ELEMENTS 

AND REORDER NODES SO THAT FIRST TWO ARE ON THE BOUNDARY. 

GEQMDR (GEOMFL) 

COMPUTE ENTRIES FOR B112 AND B113 MATRICES. 

GEOMFL (LINKl (3) ) 

GENERATE THE UNIQUE ELEMENT MATRICES AND VECTORS. 

SET UP THICKNESS VECTOR ITK. 

GENERATE LENGTH * THICKNESS ARRAY IX1P2. 

GENERATE AREA * THICKNESS ARRAY I AREA. 

GENERATE B112 MATRIX. 

GENERATE B113 MATRIX. 

IF KODG .GT. 0» PRINT ELEMENT NQ.» NODES OF ELEMENT AND 
COORDINATES OF NODES FOR ELEMENTS FROM 'IBOT' TO 'ITOP'. 

AFTER THE ELEMENT LOOP IS COMPLETED » PRINT THE VECTORS AND 
MATRICES THAT WERE GENERATED IN THE ELEMENT LOOP. 

HETALC (XYCRDM) 

COMPUTE ' ALC' AS THE SHORTEST SIDE OF ALL THE ELEMENTS » 

IF IT IS NOT READ IN NAME02. 


omrs. \ iw:s i3 

OF POO?. QUALITY 


INITIALIZATION ROUTINES 


INDEX (BDINPT »LINK3) 

III ASE LOCATIONS FOR VARIABLES TO BE USED IN SOLUTION. 
SET BASE ADDRESSES FOR ELEMENT MATR1CEX AND VECTORS 
USED IN ' DERVBL ' ELEMENT ASSEMBLY LOOP. 


JNCINP (TBLINP) 

DISTRIBUTE U1 * TKE AND DISS. OVER JUNCTURE REGION DOMAIN. 


-NHELS T (ENTRY POINT IN FENAME) 

READ NAMELISTS 'NAMEOl' AND 'NAME02' WITHOUT AFFECTING 
SCALARS IN FENAME THAT MAY HAVE BEEN CHANGED SINCE 
CALL TO FENAME. 


NODPCP (BDINPT) 

PLACE PC BOUNDARY DATA ALONG DOMAIN BOUNDARY. 
DISTRIBUTE PC DATA OVER ENTIRE DOMAIN USING STRF. 


OUiCiliV l'* ts 

OF POOR QUALITY 


OR IH DISCRETIZATION ROUTINES 


COLS (DSCRTZ) 

COMPUTE THE NUMBER OF COLUMNS* 'LCOLS IN THE OUTPUT DISPLAY 
AND SET UP THE FOLLOWING ARRAY » 

INCOL(J) - NO. OF NODES IN COLUMN J. 

DELADD < ADDDEL) 

ADD ENTRIES TO AN INTEGER ARRAY 'NS IDF ' AT A TIME. 

DELELM (DELNOD) 

DELETE ENTRIES IN AN INTEGER ARRAY 'NSIDE' AT A TIME. 

DELETE (DSCRTZ) 

DELETE NODES THAT ARE NOT CONNECTED TO ANY ELEMENTS. 

DELNOD (ADDDEL) 

3ET UP CALL TO DELELM AND SUPPRESS ZERO ENTRIES IN A VECTOR. 
DSCRTZ (LINK2U4)) 

SET UP SOLUTION DOMAIN DISCRETIZATION. 

SET UP SCALE FACTORS * GENERATE NODES AND ELEMENTS* COMPUTE ROW 
AND COLUMN KEYS FOR SOLUTION AND PRINTOUT, AND SCALE INPUT 
COORDINATES BY REFL. 


THE FOLLOWING KEYWORDS ARE USED TO TRANSFER FLOW - 


TYPE I - THIS TYPE OF INPUT IS FOR RECTANGULAR GEOMETRY. 


KEYWORD 

FORMAT 

COLS 

VX2SCL 

A8 

1-8 

VX1SCL 

A8 

1-8 


♦ DESCRIPTION 

READ THE FOLLOWING IN FREE FORMAT 
XO - START OF X2 GEOMETRY. 

NDIV1 - NO. OF DIVISIONS IN FIRST 
INTERVAL. 

X2 - X2 POS. AT END OF INTERVAL. 

PR1 - PROGRESSION RATIO FOR SPACING 

NDIV2 - NO. OF DIV. IN 2ND INTERVAL. 

X2 - X2 POS, AT END OF 2ND I NT. 

PR2 - PROG. RATIO FOR SPACING. 

. . . CONTINUE WITH NDIV3 X3 PR3* 

UNTIL A SCAN DELIMITER 'T' OR A BLANK 
CARD IS ENCOUNTERED. 

SAME AS VX2SCL BUT FOR DIRECTION 2. 




P*> 1 / 


• '/ . •!,’ fl» 

OF POOR QUALITY 


GRIP DISCRETIZATION ROUTINES 


DSCRTZ CQNT ♦ TYPE I INPUT* 


KEYWORD 


NDECRD 

N1 


FORMAT 


A8 

FREE 


COLS 


1-8 


descriptJon 


(AFTER COL. 0 ON NDECRD CARD. 

FOR N1 E'Q -It A RECTANGULAR GRID IS FORMED 
FROM THE VX1SCL AND VX2SCL INPUT. 

FOR N1 NE -If READ IN RECTANGULAR MESHES f 
• * PER CARD IN FREE FORMAT. 

TERMINATE READ WITH A BLANK CARD. 
NRLf NRUf NCLf NCU 
NRL - LOWER ROW NUMBER 
NRU - UPPER ROW NUMBER. 

NCL - LOWER COLUMN NUMBER. 

NCU - UPPER COLUMN NUMBER. 

E. G. 

3 8f 1 6 

PRODUCES A MESH OF NODES FROM ROW 3 
THRU ROW 8 AND FROM COLUMNS 1 THRU 6 
WITH SCALE FACTORS GENERATED 
ACCORDING TO THE PRUJ'S. 



HO i 

i 

1 

s N1 

1 1 

FREE 

I N2 

t 

FREE 

i N3 

[ FREE 

N4 

FREE i 

DONE 

A8 


1-8 


ON 'ELEM 


1-8 


GENERATE ELEMENTS FROM 
RECTANGULAR MESH. 


CARD +1 = ADD ELEMENTS IN ELEM. 

-1 = DELETE ELEMENTS IN ELEM. 
TURN DIAGONALS BELOW ROW N3. 

TURN DIAGONALS RIGHT OF COL. N4. 

SCALE XI COR BY XSCALE 

SCALE X2CO R BY YSCALE AND RETURN. 


* * * NOTE * * $ VX1SCL IS OF‘TIONAL FOR ID ELEMENTS. 


CREATE A GRID 41 ROWS BY 25 COLUMNS 
LINK2 14 T CALL DSCRTZ 
VX1SCL 

O.Of 12 1.0 1 *25f 12 2.0 0.8 T 
VX2SCL 

O.Of 20 2.0 l.lf 20 S.O 1.0 T 
NDECRD 

1 41 f 1 25 0 T 

ELEM 
DONE 


GRID DISCRETIZATION ROUTINES 


cr«n ? . w- * ® 

OF pfJOU QUALITY 


PSCRTZ CQNT. TYPE II 


THIS TYPE OF INPUT IS TO BE USED FOR 
NQN-RECT ANGULAR GEOMETRIES. 


KEYWORD 

FORMAT 

COLS. 

DESCRIPTION 

STYPE 

A8 

1-8 

NO. OF SIDES / SUPER ELEMENT. 

SELCN 

!* 8 

1-8 

NODE CONNECTIONS FOR SUPER ELEMENTS. 
EIGHT (8) / SET. 

N1 . N2 . N3» N4> N5» N6> N7» N8 J 
N1 AND N2 DETERMINE DIRECTION 1. 1 
N2 AND N3 DETERMINE DIRECTION 2. f 
N7 = N8 « 0 FOR TI ANGLES. j 

NETA 

A8 

00 

1 

H 

READ NO. OF DIVISIONS IN DIR. 1 | 

OF SUPER ELEMENTS. j 

NEPS 

A8 

1-8 

READ NO. OF DIVISIONS IN DIR. 2 ! 

OF SUPER ELEMENTS. j 

DEPVAR 
DEP. VAR. 

A8 

LIST 

1-8 

DATA INPUT ON SUPER NODE SEQUENCE. i 

NVAR OF THESE ARE TO BE ENTERED. ! 


FREE 

AFTER 

COL. 8 LIST OF DEP. VARIABLES 
TO BE ENTERED. 

SET STARTS WITH 89. 90 FOR ETA. EPS. 
COORDINATES. AFTER THAT (IF MORE) 

THE ORDER IS UNIMPORTANT. 


A CHECK IS MADE TO SEE IF NSNODE * NVAR EQUALS THE 
NUMBER OF ENTRIES IN THIS SET. IF THEY ARE NQTr 
THE DATA IS PRINTED OUT AND EXIT IS CALLED. 

ANY OTHER KEYWORD ENCOUNTERED WILL CAUSE A RETURN FROM DSCRTZ. 
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*-5 

OF FOG A quality 


GRID DISCRETIZAT I ON ROUTINES 
PSCRTZ CON T ♦ TYPE II CONT* 

E. G. 

REAP THE FOLLOWING IN NAMELIST 'NAME it ' 

NIZS = 250* NSNOPE = 19 * NSELEM * 2 * NVAR * 3* NODES « 600* JPR 


ADD THE FOLLOWI NG CARPS AFTER / FEOIMN ' AND 'I PINT' I NPUT PATA« 


LINK4 

9 

T 

CALL POTENT TO SET 

UP STORAGE FOR DISCRETIZER* 

PSIBD 

1 

-2 

T 

SUPER ELEMENT 

2 HAS 

DIAGONALS OPPOSITE 

OF 1. 

LINK2 

14 

T 

CALL PSCRTZ 

TO READ IN SUPER ELEMENT 

DATA* 

NETA 

20 

20 

T 

20 ELEMENTS 

ALONG 

ETA SIDE OF S*E. 1 

AND 2. 

NEPS 

12 

12 

T 

12 ELEMENTS 

ALONG 

EPS SIDE OF S.E. 1 

AND 2* 


STYPE 

4 A T S* E* 1 AND 2 ARE BOTH QUADRILATERALS * 


SELCN 

9 10 3 1 12 13 19 2* 

1 3 14 15 19 16 17 18 T 1ST 8 ARE S*E. 1 CONNECTIONS* ETC* 
PEP VAR 289 290 1248 T TRANS* COORD** NORMAL COORD** U1 VEL 

0*0 0*9 0*02 5*0*0 0*1 0*1 0*0 1.2 0*9 *02 0*0 1*101 1*2 

1*101 1*2 

0*0 0*0 *02 5*0.0 0*0 *02 3*0*0 2*0*1 4*0.0* 

2*0*0 *758 6*0.0 *783 3*0*0 .783 0*0 
4*0.0 T 1ST SET * TRANSVERSE* 2ND SET = NORMAL* 3RD SET = U1 
PONE 


r 


Ot 

OF 


& 


POOR QUALITY 


GRID DISCRETIZATION ROUTINES 


ELEM (DSCRTZ) 

GENERATE ELEMENTS AS A FUNCTION OF NODE COORDINATE INPUT. 

USED PRIMARILY FOR A RECTANGULAR DOMAIN. 

NDECRD < DSCRTZ) 

GENERATE NODE COORDINATES FROM SUPER ELEMENT INPUT DATA. 

NODELM (LINK1 (2) ) 

SET UP THE ARRAYS 'IELS' AND ' IELEM' TO STORE THE NUMBER OF 
ELEMENTS PER NODE AND A LIST OF ELEMENTS CONNECTED TO EACH NODE. 
ALSO 'AVTHK' AND 'AVAREA'* 

-PRATIO (DSCRTZ) 

COMPUTE NODES USING PROGRESSION RATIO AND END POINTS. 

REORDR (FINDBE) 

REORDER THE NODES OF AN ELEMENT SO THAT THE FIRST TWO WILL BE 
BOUNDARY NODE SPECIFICATION IN THE 'IBORD' VECTOR MUST BE 
COUNTER-CLOCKWISE . 

ROWS (DSCRTZ) 

COMPUTE THE NUMBER OF ROWS* 'KROW', IN THE OUTPUT DISPLAY 
AND SET UP THE FOLLOWING ARRAYS t 
INROU(I) - NO. OF NODES IN ROW I. 

INDRW(I) - COLUMN NUMBERS OF NODES IN ROW I. 

INDEX (J) - ROW NUMBERS OF NODES IN COLUMN J. 

NOCOL(I) - STARTING COLUMN NO. FOR ROW I. 

XYSCAL (DSCRTZ) 

COMPRESS A VECTOR OF NUMBERS 'XI' BY SCALE FACT 'SCFT'. 

FIND 'XYD « MAX (XI > - MIN(Xl) * SCFT' 

IF TWO ADJACENT POINTS OF ARRAY 'XI' ARE WITHIN 'XYD' OF EACH 
OTHER t SET THE UPPER VALUE EQUAL TO THE LOWER VALUE. 
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o»* poo-a gu/u. ity. 


II GRID DISCRETIZATION ROUTINES 


REFINED DISCRETIZATION ROUTINES. (NON-RECT ANGULAR GRID) 
(REFINE) 

GENERATE BOUNDARY CONDITION VECTOR FROM 
SURER ELEMENT BOUNDARY CONDITIONS. 

CARRAY (REFINE) 

SET UP SUPER ELEMENT COORDINATE ARRAY. 


£LKEY (REFINE) 

GENERATES REFINED GRID FINITE ELEMENT CONNECTION DATA. 

EXTRCT (REFINE) 

EXTRACTS REFINED GRID FINITE ELEMENT DATA FROM REFINED NODAL DATA. 
EXT ID (BND) 

SAME AS EXTRCT BUT FOR BOUNDARY FINITE ELEMENTS. 

MESH (DSCRTZf POTENT) 

GENERATE GRID MESH FOR VARIOUS TYPES OF INPUT. 


.POTENT (LINK4(9) ) 

S m*?r UF TE ” P0RARY STORAGE LOCATIONS IN IZ ARRAY FOR 


■QUADR (REFINE) 

PERFORMS BI-QUADRATIC TRANSFORMATION OF QUADRILATERAL 
SUBREGION DATA AND GENERATES REFINED GRID DATA. 

IF KDUMP .EQ. 1 t PRINT PRATIOS AND RSCALES. 


REFINE (MESH) 

PERFORMS GRID REFINEMENT OVER TUO-DIMENSIONAL SOLUTION DOMAIN. 
SELBND (REFINE) 

RENUMBERS THE SUBREGION GENERATED BOUNDARY GRIDPOINTS. 


SELINK (REFINE) 

FORMS SUBREGION CONNECTION TABLE. 


SETOSN (REFINE) 

FERFORMS SUBREGION ELEMENT TO SUBREGION NODE DATA TRANSFER. 
TRIANQ (REFINE) 

iKrarjjs: . 0F tr, “* 8ubreoi °" 


■XCOORD (REFINE) 

TRANSFORMS SUBREGION COORDINATES TO RECTANGULAR CARTESIAN. 


J 



INTEGRATION ROUTINES 


: t ........ V* 

UK KOO.I QUALITY 


BCONDT < LINKK7) > 

DUMMY ROUTINE FOR THE PRESENT. 

BLSF'RN ( SCHPRN ) 

PLACES CALLS TO PPRES AND STRF TO GENERATE PHI AND PP. 

XQNTES (LINK2( 4 ) ) 

RUNNING SMOOTH CONTINUITY EQUATION SOLVER TO COMPUTE U2 UP 
COLUMNS OF NODES AFTER VST ART HAS BEEN REACHED. 

IF IURIT .GT. O. PRINT DEBUG INFORMATION AT PRINT POINT. 

DERVBL (LINKK8) ) 

FORM THE DERIVATIVE OF THE ORDINARY DIFFERENTIAL EQUATION FIRST 
ON Ul-VELOCITY (GLOBAL CONTINUITY) AND OTHER DEPENDENT 
VARIABLES INCLUDING TURBULENT KINETIC ENERGY » DISSIPATION FUNCTION. 
LONGITUDINAL AND LATERAL MOMENTUM » IF DESIRED. 

IF KODS .GT. O. PRINT DEBUG INFORMATION AT PRINT POINT 
FOR ELEMENTS CONTAINING NODES 'IBOT' AND 'ITOP'. 

■DERVPX (CONTES) 

3-POINT FORWARD DIFF. FORMULA TO COMPUTE DERIVATIVE 
IN DOWNSTREAM DIRECTION. 

DFCFBL (LINK2(2) > 

COMPUTE TURBULENT VISCOSITIES FOR DEPENDENT VARIABLES FROM 
1* TKE - DISSIPATION EQUATIONS OR 

2. MIXING LENGTH THEORY. OR 

3. COMBINATION OF BOTH TKE - DISS. AND MLT . 

IF I WRIT .GT. 0. PRINT DEBUG INFORMATION AT PRINT POINT. 

DFCFNS (LINK2( 1 ) ) 

COMPUTE LAMINAR VISCOSITY ACCORDING TO TEMP. AT NODES 
USING SUTHERLAND'S LAW. 

DRHOBL <LINK2<?>> 

“CALL IF IGAS ■ 0 IN NAME01 

COMPUTES THE TEMPERATURE AND DENSITY USING A SIMPLIFIED 
ENERGY EQUATION* IS0ENERGET1C MIXING FLOW WITH 2 SPECIES* 


-OETFSL ( DFCFBL tTRBTHK) 

FIND BOUNDARY LAYER THICKNESS. DELTA. AND NODE AT WHICH IT OCCURS. 


r . f » I*.-. - f f* 
y , i J « ««•<* 

OF l-OOki QJAUlY 


INTEGRATION ROUTINES 


GET RPR (LINKK4) t NODPPR) 

COMPUTE PRESSURE GRADIENT FROM THREE POINT I NT » SCHEME USING 
INUISCID PRESSURES . 

IF IDIFRT «GT • 0» PRINT DEBUG INFORMATION AT PRINT POINT. 

JNCPPR < NODPPR ) 

GENERATE PRESSURE AND PRESSURE GRADIENT VECTORS FROM TABLE DATA. 
NODPPR ( IMPLCT ) 

SET MULTIPLIERS FOR PP RHS TERM 2. 

CALL DPSISG TO DETERMINE IF THIS IS A DATA PRINT/PUNCH POINT. 

CALL JNCPPR TO COMPUTE PRESS. AND PRESSURE GRADIANT FOR NEXT STEP. 

NUGEOM < LINKS ( 1 ) ) 

COMPUTE H21» G22r G23» FI AND F2 FOR NORMAL DIRECTION* 

H31 » H32 1 H33» G1 AND G2 FOR CROSS DIRECTION 
WHEN RUNNING VARIABLE GEOMETRY. 

PPRES ( BLSPRN ) 

COMPUTE PHI. 

COMPUTE PERTURBATION PRESSURE. 

COMPUTE GRADIENT PHI CONTRIBUTION TO U2> U3 MOMENTUM EQNS. 

IF K0D6 .GT. 0» PRINT DEBUG INFORMATION AT PRINT POINT 
FOR ELEMENTS CONTAINING NODES 'IBOT' AND 'ITOP'. 

PRSGRD (LINKl(S)) 

COMPUTE AXIAL PRESSURE GRADIENT. 

IF IDIFRT .GT. 0» PRINT DEBUG INFORMATION AT PRINT POINT. 

QKNINT (BD1NPT) 

INTEGRATION CONTROL ROUTINE TO TRANSFER CALL TO OUTPUT 
PACKAGE AT PRINT STATIONS. 

RNLDST (BLSPRN) 

COMPUTE REYNOLD STRESSES FROM TKE» DISS. FUNCTIONS. 

IF KDUMP .EQ. 1 .AND. KOUNT »LE« 2» PRINT DEBUG INFORMATION. 

SCHPRN <LINK5(S> » IMPSLVt ZZZZZZ) 

USUALLY A USER PROVIDED SUBROUTINE TO CONTROL FLOW OF PROGRAM 
AT THE END OF EACH ITERATION. 

SETDIF <LINKS(&) * IMPSLV) 

COMPUTE EFFECTIVE VISCOSITY FOR DEP. VAR. USING DFCFNS t DFCFBL. 


ORIGINAL PASS IS 

OF POOR QUALITY 


•V 


INTEGRATION ROUTINES 


SOLVER ( IMPLCT.QKNINT ) 

CONTROL ROUTINE TO DETERMINE WHICH INTEGRATION SCHEME TO USE* 

TAUW (WLFLXS) 

COMPUTE SKIN FRICTION USING PATANKER AND SPALDING 
OR LUDWIEG - TILLMAN. 

IF IDIFRT *GT* 0* PRINT DEBUG INFORMATION AT PRINT POINT* 

ICHECK ( DERVDX t GETPPR ) SUB* TCHECK ( TL, KOBE ) 

CHECK FOR ENOUGH POINTS TO COMPUTE THREE POINT INTEGRATION* 

TRAPIN ( PRSGRD * TRBTHK ) 

TRAPEZO I DAL INTEGRATION. 


TRBTHK < IMPSLV) 

COMPUTE AND PRINT INTEGRAL PARAMETERS DELTA-* AND THETA* 


WLFLXS <LINK2(3)> 

co s^ 1 ^ e ,o 8 is i 2lS^ c ^^!: ribution and heat transfer 

IF IDIFRT *GT, 0» PRINT DEBUG INFORMATION AT PRINT POINT. 


XYCRDM (LINK2( 13) > 

GENERATE VECTORS FOR GRID OUTPUT ROWS AND COLUMNS » 
®**ALE COORDINATES WHEN RUNNING VARIABLE GEOMETRY, 


ZZZZZZ (BDINPT. IMPSLV) 

^"“the^end^of^each^i terat i on?^ T IRE T0 C0NTR0L PR0BRflM 


FLOW AT 



IMPLICIT SOLVER ROUTINES 


ASMMAT < DERVBL > 

ASSEMBLE TRI-PENTA-ETC. DIAGONAL MATRIX. 

BANDET < IMPSLV) 

“liECdjHPOSE JACOBIAN INTO L-U DECOMPOSITION. 

BANDSL (IMPSLV) 

SOLVE FOR DQ USING L-U DECOMPOSITION AND BACK SUBSTITUTION. 

FILTER ( IMPSLV r SCHPRN ) SUB. FILTER < LVr KVAD t LW» KUAD ) 

AVERAGE U2 r U3 DELTA-0 IN IMPSLV. 

AVERAGE U2r U3 VELOCITIES OVER ELEMENTS. 

IMPLCT (LINK4 ( 2) ) 

IMPLICIT SOLUTION CONTROL ROUTINE. 

FOR CERTAIN ERRORS r 'TIME' IS SET TO -1.0 AND SELECTED DEBUG 
INFORMATION IS PRINTED ALONG WITH A NODE MAP AND STANDARD PRINT 

IMPSLV (DERVBL) 

IMPLICIT SOLVER ROUTINE USING BANDET AND BANDSL. 

IF NR .GE. NELPAS r PRINT DEBUG INFORMATION AT PRINT POINT; 

AND IF NIT .GT. 0» PRINT JACOBIAN AND L/U DECOMPOSITION. 

SET IMP ( IMPSLV tSTCODE) 

STORE THE J-TH LEVEL OF DEPENDENT VARIABLES FOR NEXT STEP. 
STCODE (BDINPTf DERVBL) 

“pRETOCT SOLUTION FOR CERTAIN VARIABLES (IF CALLED FOR). 


BANDED CHOLESKY SOLUTION ROUTINES. 


ASMSQ (STRF) 

BOOLEAN ASSEMBLY OF SQUARE SYMMETRIC MATRIX. 

BANCHO (STRF) 

SYMMETRIC BANDED CHOLESKY LINEAR ALGEBRAIC EQUATION SOLVER. 

STRF ( L INK2 ( 7 > t BLSPRN > 

""IMPLICIT EQUATION SOLVER USING BANDED CHOLESKY BACK SUBSTITUTION 
SOLVES FOR? PHI AND PP. 

AT INITIALIZATION - DISTRIBUTES PC OVER DOMAIN IN NODPCP. 

IF IPWRIT .GT. Or PRINT DEBUG INFORMATION FOR ELEMENTS 
FROM ' IBOT' TO 'ITOP'. 

SUMKEY (LINK3(5>> 

"•g gHPU TE ELEMENT LINKING KEYS FOR SOLVER ROUTINE BANCHO. 


ORIGINAL PAGE IS 

ai i JM ITV 


OUTPUT ROUTINES 


CALORD < BDINPT > 

CALL ORDER OF ROUTINES USED FOR VARIABLE 
PRINT INTEGRATION VARIABLE NOS. AND DEPENDENT VARIABLE^NOS 
PRINT LIST OF PARAMETERS TO BE PRINTED IN OUTPUT ROUTINE? 
PRINT PLOT INFORMATION* KUUTlNfc. 


CQMQC (BDINPT) 

PRINT THE COMOC SYMBOL ON TWO PAGES ALONG WITH ASSOC. TITLE CARDS. 
DPSISQ < NODPPR ) 

DETERMINE IF THIS IS A PRINT/PUNCH DATA STATION AND SET 
APPROPRIATE FLAGS. 


DRVBUG (DERVBL) 

DEBUG PRINT ROUTINE FOR DERVBL. 

THIS ROUTINE IS CALLED AT A PRINT POINT WHEN 

K0D5 .GT. 0 AND AN ELEMENT CONTAINS NODE 'IBOT' OR 'ITOP' 


-FEPLOT < STOUT 1) 

GENERATE DATA TO BE USED FOR PLOT PACKAGES. 

TH IL.?° U Ii NE IS CALLED AT THE END OR A PRINT POINT WHEN: 
KPNT .GT. 0 AND KOUT .GT. 0 AND KPLVAR .GT. 0 


ICQND (BDINPT) 

PRINT INTEGER AND REAL SCALAR CONDITIONS. 
I ARRAY ( 1 ) - lARRAY(SOO) 

RARRAY ( 1 ) - RARRAY(SOO) 


fl-UTNOP SUB. OUTNOD < NN» IRRAYr TITLE > 

PRINT AN INTEGER ARRAY ALONG WITH A 32 CHARACTER TITLE. 


OUTPG (GEOMFL) 

PRINT THE ELEMENT NO. AND NODE CONNECTIONS AND NODE 
COORDINATES FROM THE GEOMETRY ROUTINE 'GEOMFL' ♦ 

THIS ROUTINE IS CALLED WHEN KODG .GT. 0 AND THE ELEMENT 
NO. IS IN THE RANGE FROM 'IBOT' TO 'ITOP'. 


Oi t QUAUrV 


OUTPUT ROUTINES 


.OUTVEC SUB. OUTVEC ( NN» ARRAY » TITLE > 

PRINT A REAL ARRAY ALONG WITH A 32 CHARACTER TITLE. 

PBLANK (REOUTP) SUB. PBLANK ( N» P ) 

INSERT 'N' BLANKS IN THE OUTPUT VECTOR 'P' . 

-PLILNK (REOUTP) 

CONVERT A FLOATING POINT NUMBER INTO 'A' FORMAT. 

PRINTA SUB. PRINT A < P. NPLC. VAR ) 

PRINT A LIST OF 'NPLC' REAL NUMBERS ' VAR ' IN 'A' FORMAT. 

REOUTP <LINK2<5>> 

PRINT THE ARRAY GEOMETRY AND NODE NUMBERS IN A PATTERN THAT 
RESEMBLES PROBLEM GEOMETRY. 

(FEOUTP ) <LINK2( 6) ) FEOUTP IS A TRANSFER POINT IN REOUTP 

WHEN KROW .LT* 0. 

PRINT OUTPUT PARAMETERS IN A PATTERN THAT RESEMBLES 
PROBLEM GEOMETRY. 

IF MAX. SCALE FACTOR EXCEEDS 'NSM' (DEF. = 10). 

TERMINATE THE PROGRAM. 

IF OUTPUT PRINT NO. ' KOUNT ' . EXCEEDS PRINT LIMIT 'LPRINT' 
TERMINATE THE PROBLEM. 

IF NEWPRT .GT. 0. PRINT SCALAR OUTPUT ' NEWPRT ' COLUMNS ACROSS. 
IF NMOUT .EU. 2. PRINT VECTOR OUTPUT WITH A CALL TO OUTVEC. 

RITE (LINK3(2) ) SUB. RITE ( NMB. KEY. TITLE. NTITL ) 

"CORPUTE ' NUMBER * <KEY-1>»10 + NMB '. 

GO TO STATEMENT ACCORDING TO VALUE OF 'NUMBER'. 

IF 'NUMBER' IS OUT OF RANGE. WRITE TITLE INFORMATION. 

SCALEV (REOUTP) 

CALL SCALE ROUTINE FOR UP TO 10 OUTPUT VARIABLES AT A TIME. 

SETSCL (SCALEV) SUB. SETSCL ( ARRAY. ILIM, IT ) 

SET SCALE FACTOR FOR AN ARRAY OF REAL NUMBERS 
AND NORMALIZE THE ARRAY BY 10**IT. 

STOUT 1 (REOUTP) SUB. STOUT 1 ( MTM. NVOUT ) 

DIMENSIONALIZE OUTPUT VARIABLES FOR DISPLAY PURPOSES. 


.#*•'« M page is 


DUMMY ROUTINES 


GET ADD (FEDIMN) SUB* GETADD < 1VAR'» IADD > 

INSTALLATION ROUTINE TO GET MACHINE ADDRESS 
OF VARIABLE 'IVAR' AND STORE IT IN 'IADD'* 

GETDAT <RDATE»RITE) SUB* GETDAT ( DATE > 

INSTALLATION ROUTINE TO GET CURRENT DATE ON SYSTEM 
AND STORE IT IN 'DATE 7 ♦ < DIMENSION OF DATE IS 2 ). 

6ETTIM 

INSTALLATION ROUTINE TO GET CURRENT CPU TIME. 

L1NK3A (LINK3< A) > 

USER ROUTINE TO BE USED AS NEEDED* 

PDUMP (BDINPTf FEDIMN) SUB. PDUMP < A» B* N > 

INSTALLATION ROUTINE TO DUMP CORE FROM 'A' TO 'B' 

IN FORMAT 'N'. 

RDATE 

PLACES CALL TO GETDAT. 


TIMETK SUB* TIMETK ( T > 

INSTALLATION ROUTINE TO STORE CURRENT CPU TIME INTO 'T 


Sirs* +**^1 a'’ t ? 

(mt.'rfikviW tv,v.« 

OF POOR QUALITY 


V 


THE FOLLOWING CONTAINS A SEQUENTIAL LIST OF THE ENTRIES 
IN ' I ARRAY' ♦ THE REFERENCES ARE TO THE SUBROUTINES IN 
WHICH THEY ARE USED ♦ IF NO REFERENCE IS GIVEN* THEN THE 
VARIABLE IS USED IN NUMEROUS ROUTINES. 


SEQUENTIAL IARRAY ENTRIES 


NO. 

ENTRY 

REFERENCES 

1 

ND 

DFCFNS* 

6 

KODG 

GEOMFL » 

7 

K0D5 

BCONDT * DERVBL » INDEX* PPRES* STCODE* 

8 

KPRINT 

CONTES* 

14 

NELEM 


16 

NNODE 

i 

19 

NOE 

♦ 

20 

NPRNT 

LINK3* 

22 

NC 

LINK2* 

23 

NB 

LINK2* 

26 

KOUNT 

BDINPT * LINK2* QKNINT * RNLDST * 

27 

NSKIP 

BDINPT* BCONDT* FEDIMN* GETBCD » POTENT* > 

28 

IPASS 

LINK2* DERVBL* 1MPLCT * STCODE* STRF * 

30 

NP 

! 

31 

NEQ 

\ 

1 

32 

IWLSEP 

XYCRDM* l 

34 

LPRINT 

LINK2* QKNINT* 1 

37 

NEWPRT 

REOUTP * j 

38 

NOUTVC 

REOUTP » 

40 

NOUTS 

RFOUTP* 

42 

INPUT 

E.INPT* DESCRP* GETBND * READV1 * 

43 

NEQADD 

BLS1US* CONTES* DFCFBL* IMPLCT » IMPSLV* 

46 

NF 

LINK2* DESCRP* FEDIMN* 

47 

LG 


50 

LCOL 

i 

» 

52 

KROW 

S3 

NDB 

LINK2t IMPLCT* ! 

54 

NHHALF 

LINK2p IMPLCT* ; 

55 

NODE 

1 

58 

NEGKNN 

♦ 

i 

59 

NCPTAB 

CPINIT* FEDIMN* i 

60 

NMBOUT 

LXNK2* CALORD* DESCRP* FEDIMN* ST0UT1* i 

61 

KDUMP 

BDINPT* LINK2* LINK3* COMOC* DESCRP* DIMEN* FEDIMN* ] 

62 

NTITL 

64 

NSM 

LINK2* 

68 

NI 

LINK3* NBNDRY * « 

86 

KPNT 

- 


35 



1<r / 

ORIGINAL PAGE IS 

OF POOR QUALITY 


SEQUENTIAL I ARRAY ENTRIES 


NO. ENTRY REFERENCES 

88 NLINE COHOC » DFCFBL* TRBTHK * 

90 NYY LINKS* LINK3* DERVBL » DSCRTZ* FEDIMN* INDEX* 

JNCPPR * NODPCP* PPRES* STCQDE* STOUT1* STRF » 

91 NZZ FEDIMN* INDEX* STOUT1* 

92 IZSIZE MAIN* FEDIMN* POTENT* RITE* STRF, 

94 NIZS FEDIMN, IMPLCT* 

96 IBASE BDINF'T » LINKS* DSCRTZ* STOUT1* 

97 ITKE DFCFBL* 

99 IBL BDINPT » LINK1* LINKS* L1NK3, BCONDT, COMOC, 

DFCFNS* DIMEN, FEDIMN* FILTER* GEOMFL* IMPLCT* 
IMPSLV, RITE* STRF* 

100 IREND BDINPT, LINK2* LINKS* FEDIMN* POTENT* RITE* STRF, 
102 KOUT BDINPT* LINK2* CALORD, IMPLCT* QKNINT* ST0UT1* 

105 IPTSPL TAUW* 

107 NE1E2 DFCFBL* JNCINP* 

108 NCNADD FILTER* IMPLCT* PPRES* 

113 KPLVAR BDINPT* LINKS* CALORD* FEDIMN* ST0UT1* 

122 I WRIT CONTES* DFCFBL* 6ETFSL* IMPLCT* TAUW* WLFLXS, 

125 NCALLS LINKS* CALORD* FEDIMN* 

127 IDIFRT LINK1* GETF’PR* IMPLCT, JNCPPR* PRSGRD* 

131 NBORD BDINPT* 

132 IPURIT LINK3* STRF* 

133 I TOP BCONDT* DERVBL * GEOMFL* IMPLCT* PPRES* STRF, 

142 NOUTPR LINKS* FEDIMN* 

151 KWFLXS WLFLXS* 

152 INITKE DFCFBL, 

153 NFUNCH DIMEN* DSCRTZ* FINDBE* 

157 NPASS IMPLCT* QKNINT* 

161 NPVSX BDINPT* BRDSHU* DIMEN* DPSISQ* FEDIMN* GETPPR* 

JNCINP* JNCPPR* NODPCP* PBLANK* 

162 IEPSET DFCFBL* 

166 N3DPNS BDINPT* LINK!* COMOC* CONTES* DERVBL* IMPLCT* 

RITE* STCODE* 

167 KNTPAS DPSISQ* IMPLCT* 

169 KCDC BDINPT* 

173 NCOMOC COMOC* DESCRP* RITE* 

174 NCOMTD BDINPT* COMOC* DESCRP* RITE, 

177 NU2P0S FEDIMN* NWOEOM* 

178 NU3P0S CONTES, FEDIMN, NWOEOM* 

179 LOC DIMEN* GETFSL, GETPPR, LOOK* LOOKAV* 

186 ITDA CALORD* TRBTHK* 

187 ITDB ST0UT1 * 




OF POOR QUALi i ( 


SEQUENT I AL I ARRAY ENTRIES 


NO. 

ENTRY 

REFERENCES 

189 

NSTD 

BDINPT » FEDIMN# 

190 

NMOUT 

L1NK2# LOOK# 

191 

NM 


192 

N2M 

DERVBL# FEDIMN# GEOMFL# RNLDST » STRF » 

193 

NH2 

DERVBLf DRVBUQ t FEDIMNt GEOMFL t PPRESt j 

RNLDST » STRF# 

194 

NDP 

BDINPT# LINK3# FEDINN# IMPLCT# INDEX# STCODE# 

195 

KFXBND 

WLFLXS » 

196 

1TUALL 

ULFLXS# 

197 

NTPRNT 

LINK2# TRBTHK# 

203 

NTCRDM 

POTENT# XYCRDH# 

205 

ITIMER 

LINK1# LINK2# LINK3# LINK4# LINKS# SOLVER# 

206 

NMDL 

DERVBL# DRVBUG# FEDIMN# RNLDST# STCODE# 

208 

NVRHS 

LINK2# BCONDT# NODPCP# PPRES » STRF# 

210 

NC0HP6 

BDINPT# LINK4# FEDIMN# PRSGRD # 

212 

MLTRHS 

FEDIMN# PPRES# STRF# * 

215 

ISUPRS 

BDINPT# LINK2# j 

225 

NVRH 

LINK2# BCONDT# NODPCP# PPRES# STRF# 1 

233 

NODES 

POTENT # j 

250 

IPHI 

BDINPT# LINK3# GEOMFL# RITE# ' 

251 

NSNODE 

DSCRTZ » MESH# POTENT# \ 

252 

NSELEM 

DSCRTZ# MESH# POTENT# 

260 

NVAR 

DSCRTZ# MESH# POTENT# * 

286 

NEQAV2 

IMPLCT# 

287 

NEQAV3 

IMPLCT# 

291 

JPR 

QUADR# 

295 

I PLOT V 

CALORD# ! 

301 

NIMPLT 


303 

NITER 


304 

NELPAS 

i 

305 

NBAND 

ASMMAT# DERVBL# IMPSLV# INDEX# 5 

306 

NH8JAC 

DERVBL# IMPSLV# INDEX# l 

307 

NCONV 

BLSZZZ# DERVBL# DFCFBL » IMPSLV# PPRES# 1 

PRSGRD# TRBTHK# • 

310 

NFMROU 

DERVBL# IMPSLV# INDEX# STCODE# ; 

312 

IUONLY 

DPS1SQ# IMPSLV# PPRES# 1 

313 

NR 

IMPSLV# JNCPPR# STCODE# ! 

316 

NIT 

IMPSLV# ; 

317 

NDBGPT 

DERVBL# IMPLCT# IMPSLV# PPRES# STCODE# STRF# 

320 

ICHI 

IMPLCT# ; 

322 

I START 

IMPLCT# 

323 

NLAST 

IMPLCT# 

f 



QfiteittM. i ' 

OF POOR Q‘ 


SEQUENTIAL I ARRAY ENTRIES 


ENTRY REFERENCES 


BLSPRN 
BLSXUS 
IMPSLV 
DERVBL 
PRSGRD 
JNCINF* 
IMPSLV 
MATSUM 
HATSUH 
DPS ISO 
DPSISQ 
DPS1SQ 
DPSISQ 
DPSISQ 
DPSISQ 
NODPCP 
PPRES » 
BCONDT 
DPSISQ 
! JNCPPR 


FILTER* IHPLCT * IMPSLV* PPRES 
GETPPR » IMPLCT t IMPSLV* 


NODPCP »STRTCH* 


SET UAL > 
SETVAL* 


REOUTP t 


GETPPR * IMPLCT* PPRES * 
JNCPPR* NODPCP * 







cQcv’a 

OF POOR QUALITY 


t« E <6ADcS!! NG CONTAINS A SEQUENTIAL LIST OF THE ENTRIES 
IN RARRAY ♦ THE REFERENCES ARE TO SUBROUTINES IN 
WHICH THEY ARf USCD. IF NO REFERENCE IS GIVEN* THEN THE 
VARIABLE IS USED IN NUMEROUS ROUTINES# 


NO* 


1 

2 

3 

4 
7 

9 

10 

13 

14 

15 


17 

21 

22 

23 

24 

26 

27 

28 

30 

31 

32 

35 

36 
38 
3? 
40 

42 

43 


ENTRY 


FACT 

ONE 

ALC 

THK 

HSINIT 
PINF 
RHOINF 
DELP 
EPS 
H 


SEQUENTIAL RARRAY ENTRIES 
REFERENCES 


LJ.NK2* BCONDT » 
TRBTHK, 


DIMEN* FEPLOT * GETPPR, JNCPPRi 


16 HMAX 


* HMIN 
RE 

TF 

TIME 

TO 

TWOPI 

UINF 

RUNIV 
CPOINF 
6 

RTC0N1 
TD 
PEDDlMl 
XMUINFj 
PED6E 
TRATIol 
TKEDGE] 
REFL 


DIMEN, 

DIMEN* 

BRDSHU 

DFCFBL 

DIMEN* 

DIMEN* 

LINK2* 

INDEX* 

DIMEN* 

LINK2* 

BCONDT 

INDEX* 

LINK2* 

QKNINT 


GEOMFL* 

* DIMEN* 

* DFCFNS 
IMPLCT * 
IMPLCT * 
BLSPRN* 
JNCPPR, 
IMPLCT* 
BLSIUS* 

* BLSIUS 
SETDIF * 
CONTES* 


GETPPR, 
DRHOBL , 


JNCPPR* PBLANK, 
WLFLXS* 


* PPRES* STRF * 

* DIMEN, DPSISQ, IMPLCT, IMPSLV. 


DRHOBL 

* DIMEN 
QKNINT 
IMPSLV 
DERVDX* 

LOOKAV* 

IMPSLV* 

DIMEN* IMPLCT* PRSGRD* 

* CONTES* DERVBL * DERVNS* DFCFBL* DIMEN 
TAUW* TRBTHK* 

DERVBL* DIMEN* DPSISQ* IMPLCT* IMPSLV* 


BLSIUS* CONTES* 
NWGEOM* STCODE* 
DIMEN* GEOMFL* 
BCONDT* BLSIUS* 
JNCPPR* PRSGRD* 


DIMEN* DPSISQ* DRHOBL* IMPLCT* JNCPPR 


DIMEN* DRHOBL* 
CPINIT* DIMEN* 
BCONDT, DIMEN* 
CPINIT* DIMEN, 
DIMEN* DPSISQ, 
DIMEN * OETPPR * 
BCONDT* DFCFBL* 


CPINIT* 

WLFLXS, 


DFCFBL, DIMEN* DRHOBL* INDEX i 


PRSGRD* 

FEDIMN, 

FEDIMN* 


JNCPPR, 

DIMEN* 


DIMEN* DRHOBL* GETPPR, 
DIMEN* DRHOBL* 

DFCFBL* 

LINK2* DFCFBL* 


PRSGRD. 


PRSGRD* 

SETDIF, 

JNCPPR* 


PRSGRD. 


DIMEN* DSCRTZ* JNCPPR* PBLANK* PRATXOi 




r 


OKl&IHAL W&& & 

OF POCH QUALITY 


SEQUENTIAL R ARRAY ENTRIES 


NO. 


ENTRY 



47 

REFLRE 


48 

TIMESV 

= J 

SO 

SSINIT 


56 

RTC0N5 

■ 

58 

TOFINF 

* 

5? 

FACTMU 


60 

GAMMAF 

t • 

61 

XMACHO 


62 

CONV 

. 

63 

UEDGE 

— : 

66 

XMF 


67 

PR 


68 

EP4MD 

■j; 

70 

CONRHO 

*4 

71 

STLDVR 

y> 

Sif. 

72 

STLDTR 

73 

STLDCR 


74 

STLDEX 

- 

75 

STLCON 

V 

79 

FACTF* 

t 

80 

FACTH 

Y ‘ 

82 

PI BAR 


83 

COMPX 

; V.’ 

84 

COMPY 

' i ; 

86 

AMD 

: \ 

89 

EPS INF 

■ ", 

90 

TKEINF 


95 1 

EPTEST 

/ 

98 

XPRIME 


99 

PPRCON 

• *5 

100 

PFRIME 

■"i 
' ’j 

102 

VST ART 

! • •) 

: A 

: 1 
T 

_* T 

103 

DEPLT 

104 

VELCST 

109 

XMA 


110 

XMH 

i. 

111 

TWELVE 

! i: 

115 

AR8CAL 


124 

CON 


125 

XLAM 


REFERENCES 


DIHEN » 

IMPLCT* IMPSLV* PRSGRD » 

DIMEN* IMPLCTf QKNINT » 

CPINITf DIMENt 

CPINITt DFCFNSf DIMEN f DRHODL * FEDIMN* PRSGRDf 
DFCFBLf DIMENf ULFLXSf 
DIMEN* DRHODL* PRSGRD* 

DIMEN* DRHODL* 

BRDSHU* DIMEN* DRHODL* GETPPR* PPRES* TRDTHK* 
DIMEN* 

FEDIMN* 

PRSGRD* 

DIMEN* STRF * 

DIMEN* DFCFNS* 

DIMEN* DFCFNS* 

DIMEN* DFCFNS* 

DIMEN* DFCFNS* 

DIMEN* DFCFNS* 

BDINPT* DIMEN. STRF. 

DIMEN* 

JNCINP. 

DSCRTZ * 

DSCRTZ* 

DFCFDL » DIMEN* 

DIMEN* 

DFCFDL* DIMEN* 

DERVBL* DFCFDL* 

DIMEN* GETPPR* JNCPPR* PRSGRD* 

DIMEN* GETPPR* JNCPPR* PRSGRD* 

DIMEN* GETPPR* JNCPPR* PRSGRD* 

DLSIUS* CONTES* DIMEN* 

DIMEN* QKNINT* 

DIMEN* DIMEN* DRHODL* 

DIMEN* DRHODL* PRSGRD* 

DIMEN* 

BRDSHU* 

PRSGRD* 

DFCFDL* 

DFCFDL* GETFSL* 


\ 


I 




40 


vrr^. b 

OF POOR QUALITY 


SEQUENTIAL ft ARRAY ENTRIES 


,NO. 

entry 

12? 

SCT 

13S 

ENERGY 

136 

OSMAX 

137 

AVDP 

13d 

RUESQ 

13? 

H31 

140 

! 632 

141 

G33 

1*2 

G1 

143 

C4EDSW 

145 

E1E2SW 

154 

XMACHS 

155 

1 TSINF 

156 

AINF 

160 

CPINF 

171 

EULER 

17? 

RADCON I 

181 

C1K0RE 1 

182 

C1D0RF I 

183 

i C2D0RF 

184 

C2K0RE [ 

185 

PRDIS I] 

186 

H21 | 

187 

622 1] 

188 

623 || 

18? 

FI t 

1?0 

SLOPE ( 

l?3 

XNWGEO C 

1?6 

THETAP C 

1?8 

yplus I 

1?? 

RNULOC J 

201 

XTC T 

202 1 

UWALL T 

203 1 

JEDGEN | D 

204 1 

ENER Id 

205 1 

BELSTR Id 

206 1 

rHETA ID 

207 | 

:onvrg II 

20? | 

M.TH 1 01 


REFERENCES 


DFCFSLf FEDIMN, 

DERVBL, 

BCONDT, DFCFBL, NUGEOMi 
DFCFBL, 

TAUW» WLFLXS, 

BCONDT » DERVBL , DERVNS, 
RNLDST , STRF , XYCRDM, 
BCONDT , DERVBL, NWGEOM, 
BCONDT, DERVBL , NWGEOM, 
NWGEOM, XYCRDM, 

DFCFBL, DIMEN, 

DFCFBL, DIMEN, 

LINK3, DERVBL, 

CPIN1T, DIMEN, 

DIMEN, DRHOBL , 

CPINIT, DIMEN, 


PRSGRD, XYCRDM i 
NWGEOM, PPRES* 


PPRESi 

PPRES, 


XYCRDM, 

XYCRDM, 


DIMEN, DRHOBL, PRSGRD, 


DERVBL , 

NODPCP , 

BCONDT , 

INDEX, 

INDEX, 

INDEX, 

INDEX, 


DIMEN* 

PPRES, 

DIMEN, 


DRHOBL 

PRSGRD 

INDEX, 


GETPPR, JNCPPRi 


f BCONDT, CONTES, 
IM, PPRES, PRSGRD, 
IT, CONTES, DERVBL , 
'T, CONTES, DERVBL , 
* NWGEOM, XYCRDM, 
R> JNCPPR , LOOKAV, 
S, DFCFBL, DPS1SQ, 
R» IMPLCT, IMFSLV, 


DERVBL , 
RNLDST , 
NWGEOM i 
NWGEOM, 

NWGEOM, 
DRVBUG , 
INDEX, 


DERVNS , IMPLCT, 

STRF, XYCRDM, 

PPRES, XYCRDM, 

PPRES, PRSGRD, XYCRDM i 

NWGEOM, 


l, 

if 

> 9 

> f 


TRBTHK , 
DERVNS, 

TRBTHK, 

IMPSLV, 

TRBTHK, 







original SW- 
OP POOH QUA 


..iY 




SEQUENTIAL RARRAY ENTRIES 


NO. 

ENTRY 

REFERENCES 

346 

VLDHLT 

DERVBL* 

347 

TSCALE 

GETPPR* IHPLCT * JNCPPR* PRSGRD » 

348 

VCHULT 

DERVBL* 

3S1 

PST AG 

DIHEN * GETPPR* JNCPPR* PBLANK* PRSGRD * 

353 

GUMULT 

DERVBLt ! 

354 

UtNFX 

DIHEN* 

356 

CE 

SETDIF* 

357 

CU 

SETDIF * 

358 

EFMULT 

DFCFBL* 

359 

RETHM 

DERVBL * 

360 

RHSCAL 

RNLDST * 

361 

0S6 

FEDIHN* 

362 

0S12 

DERVBL* DERVNS* FEDIHN* GEOHFL » IHPLCT* INDEX* PPRES* 

363 

0S60 

FEDIHN* 

364 

CK 

DFCFBL* 

365 

CD 

GETFSL* INDEX* » 

368 

ESCF 

DFCFBL* I 

371 

UBAR 

DERVBL* DRHOBL* PRSGRD* j 

372 

TBAR 

DERVBL* DRHOBL* PRSGRD* ; 

373 

XMDOTC 

DERVBL* PRSGRD* ' 

374 

TAREA 

DERVBL* PRSGRD* j 

375 

0SH1SQ 

DERVBL* IHPLCT* NUGEOH * XYCRDH* 

383 

USHAX 

BCONDT * NUGEOH* XYCRDH* 1 

385 

UCMULT 

DERVBL* ! 

386 

H21L 

NUGEOH* XYCRDH* 

387 

H31L 

NUGEOH* XYCRDH* 

390 

DYNPRS 

NODPCP* 

393 

AN8 

PPRES* 

394 

T2F1X 

PPRES* 

396 

OSG 

PPRES* * 

397 

T3FIX 

PPRES* i 

398 

T2PF1X 

PPRES* 5 

399 

TU1U2P 

DERVBL* ! 

458 

U1MIN 

SETUP* ; 


« 


i 

\ 
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ORIGINAL PAG* IS 

OF POOR QUALITY 


INTEGER SCALAR VARIABLES < I ARRAY) 
CROSS REFERENCE LIST 


* * * NOTE * * * 

VALUES IN PARENTHESIS ARE DEFAULT CONDITIONS. 


ALPH ANUMER IC 1ARRAY ENTRIES 


I ARRAn 
ENTRb 


NAHE 


DESCRIPTION 


96 

261 

99 

344 

207 

327 

127 


I BASE 

IBC 

IBL 

IBLAS 

IBOT 

IDDXST 

IDIFRT 


162 IEPSET 
347 INCLFT 


(200) BASE NO. FOR IZ ENTRIES. 

NUMBER OF BOUNDARY CONDITION TYPES. 

1 = BOUNDARY LAYER PROGRAM 

EPS EXPONENT FOR U2» U3 CONVERGENCE CRITERION. 

1ST ELEMENT AT WHICH TO PRINT DEBUG INFORMATION. 

( 1 ) 

NO. OF TIMES TO PRINT INTER. OUTPUT IN LINK1 * GETPPR» 
JNCPPR AND PRSGRD. 

1 = FORCE COMPUTATION OF DISSIPATION FUNCTION. 

OUTPUT VECTOR INCREMENT FOR SETVAL AND MATSUM. 


348 


INCRGT 


158 

152 

295 

42 

28 

105 


INITCN 

INITKE 

I PLOT V 

INPUT 

IPASS 

IPTSPL 


132 

294 

100 

259 

322 


IPWRIT 

IRAT 

IREND 

ISIDE 

ISTART 


215 


1SUPRS 


1ST INPUT VECTOR INCREMENT FOR MATSUM. 

INPUT VECTOR INCREMENT FOR SETVAL. 

2ND INPUT VECTOR INCREMENT FOR MATSUM. 

INITIALIZER IN CONTES. 

1 = TKEt DISSIPATION ARE ALREADY INITIALIZED. 

DATA SET NUMBER ON WHICH TO STORE DATA FOR PLOTTING. 
(5) INPUT LOGICAL UNIT NUMBER. 

NO. OF CALLS TO DERVBL • 

(0) » USE LUDWIEG - TILLMAN FORMULA FOR TAU WALL. 

1 = USE PAT ANKER AND SPALDING'S FORMULA FOR TAU WALL. 
DEBUG CODE IN LINK3 AND STRF FOR INTERMEDIATE OUTPUT. 
SET « TO 1 WHEN N3DPNS IS FIRST TURNED ON. 

END POSITION IN 'IZ' ARRAY. 

NUMBER OF SIDES / SUPER ELEMENT. 

STARTUP CODE IN IMPLCT. 

1 « SUPPRESS PRINTOUT OFJ 
A. STARTUP IN BDINPT « 

B* NODE MAP. 

C. OUTPUT STATIONS. 


44 


Vi 


of pooh QUMUT 


ALPHANUMERIC I ARRAY ENTRIES 


I ARRAt 


NAME 


DESCRIPTION 


ENTRY! 


186 

187 

205 

97 


133 

196 

312 

122 

92 

291 


ITDA 

ITDB 

ITIMER 

ITKE 

I TOP 

ITWALL 

IUONLY 

IWRIT 

I2SIZE 

JPR 


UNIT NO. ON WHICH TO STORE INTEGRAL PARAMETER DATA 
FOR PLOTTING. 

UNIT NO. ON WHICH TO STORE 'PLOTS' DATA FOR PLOTTING. 
SET = NO. OF TIMES TO CALL TIMETK SUBROUTINE. 

NOTE 4 . TIMETK NEEDS TO BE DEFINED FOR YOUR INSTALLATION 

0 * DO NOT INTEGRATE TKE - DISS. EQUATIONS. 

1 = USE TKE - DISS. TO COMPUTE TURBULENT VISCOSITY ♦ 

2ND ELEMENT AT WHICH TO PRINT DEBUG INFORMATION. 

1 * USE DUDY FOR TAU WALL. 

(2) NO. OF STEPS UNTIL CONVERGENCE ON U1 ONLY. 

DEBUG PRINT FLAG IN CONTES , DFCFBL. GETFSL AND TAUW. 
MAXIMUM DIMENSION OF IZ VECTOR. 

1 = MIDSIDE NODES FOR ETA DIRECTION ARE 


♦ 


169 

61 

167 

6 

7 

26 

113 

86 

8 
52 

151 

50 

47 


KCDC 

KDUMP 

KNTPAS 

KODG 

K0D5 

KOUNT 

KPLVAR 

KPNT 

KPRINT 

KROW 

KWFLXS 

LCOL 

LG 


330 

214 


LMDJAC - 
LMLT 


179 

172 

376 

353 


LOC 

LOUD 

LPHI 

LPLOT 


PROGRESSION RATIOS. 

1 * RESET NLINE TO 50 AND DUMP KOBE TO 2. 

PRINT INPUT CARDS AND DATA GENERATED IN BDINPT. 

(99) MAXIMUM NO. OF INTEGRATION STEPS BETWEEN PRINTS. 
PRINT GEOMETRY OUTPUT IF .NE. 0. 

PRINT INTERMEDIATE DERVBL OUTPUT K0D5 TIMES. 

RUNNING COUNT OF OUTPUT. ( LIMITED BY LPRINT.) 

(10) NO. OF VARIABLES TO BE PLOTTED OR PUNCHED. 

1 = PRINT STATION ( SET DURING EXECUTION.) 

PRESENT VALUE OF PRINT COUNTER. 

(100) NO. OF ROWS IN DISCRETIZATION. 

(LCOL) NO. OF SLICES AT WHICH TO COMPUTE TAU WALL. 

(20) NO. OF COLUMNS IN DISCRETIZATION. 

NO. f COLS. IN SOLUTION FIELD. 

IF .NE* 0 ON INPUT » THEN CNTPTS AND CNTNDS ARE 
TO BE READ IN. 

(1) COMPUTE JACOBIAN EACH ITERATION. 

(LG) NO. OF CONTOURS FOR WHICH TO COMPUTE MIX. LENGTH 
TURBULENT VISCOSITY. 

INTERVAL NO. FOUND IN LOOK SUBROUTINE. 

(2) USE LAMINAR VISC. BELOW LOUD AND MLT FROM LOUD ON. 
+/- 1 TO EVALUATE PHI OR PP RHS IN PPRES. 

FLAG TO INITIATE PRINT OF DATA AT PLOT STATION. 


354 

352 

34 

351 

350 


LPLTPR 

LPPNCH 

LPRINT 

LPSUP 

LPUNIT 


(SEE DPSISQ AND REOUTP) 

FLAG TO ACTUATE PRINT OF DATA AT PLOT STATION (REOUTP) 
NO. OF PASSES BEFORE LPUNIT TAKES AFFECT FOR K* EPS. 
(100) LIMIT ON OUTPUT COUNT. 

NO, OF DATA STATIONS BEFORE LPLTPR TAKES AFFECT. 

UNIT NO. ON WHICH TO STORE PP* U1 r U2» U3t K> EPS DATA. 


- ... ' <•*“<* r \ r r T ♦ f r V 

r *, t *** 1 

01* S'GOK QUALITY 


ALPHANUMERIC I ARRAY ENTRIES 


X ARRAY name 

ENTR1 


DESCRIPTION 


212 MLTRHS 
377 HONE 
23 NB 
305 NBAND 
170 NBC 
131 NBORD 
22 NC 
125 NCALLS 
108 NCNADD 

173 NCOMOC 
210 NCOMPG 

174 NCOMTD 
307 NCONV 

13 NCOORD 
59 NCPTAB 

1 ND 
317 NDBGPT 
124 NDERIV 
194 NDP 

14 NELEM 
304 NELPAS 

31 NEQ 
43 NEQADD 


286 NEQAV2 

287 NEQAV3 
58 NEQKNN 
37 NEWPRT 

107 NE1E2 


46 NF 
54 NHHALF 
53 NDB 
68 NI 
301 NIMPLT 


(1) NUMBER OF RIGHT HAND SIDES TO SOLVE FOR IN STRF 
-1 IN PPRES TO ALTERNATE SION OF LF'Hl. 

(4) NO. OF CHAR* IN EACH WORD OF OUTPUT VAR. TITLE 
< 24LC0L+3) MAXIMUM BANDWIDTH OF JACOBIAN MATRIX 
MAX* NO* OF BOUNDARY COND. FOR ANY ONE DEP V VAR* 

NO* OF NODES AROUND BORDER OF DISCRETIZATION^ 

<8> NO. OF CHARACTERS IN OUTPUT FORMAT* 

° F R0UTlNES TO CALL AT END OF INTEGRATION STEP 
5n 6I JL U £lr.r? INTEGRATION AFTER NCNADD INTEGRATION STEPS 
NO* OF CARDS READ IN FOR COMOC TITLE PAGE* 

«Q 5) nr 0 ^pL E ^; E ?., T Lr R1NT FR0M PROt,UM COMMON BLOCK. 
NO* OF CARDS READ IN FOR TITLE INFORMATION 

NON-CONVERGENCE FLAG SET IN IMPSLV* 

FLAG FOR GENERATING AXI -SYMMETRIC DATA IN GEOMFi 

(1) NO* OF ENTRIES IN SPECIFIC HEAT TABLE 6E ° MFL * 
INITIALIZATION PARAMETER IN DFCFNS. 

(2) - F call R derobl. PRINT IN wectors for debugging. 

(10) SPACE ALLOCATION IN IIPINT VECTOR* 

NUMBER OF ELEMENTS IN SOLUTION* 

INTEGRATION STEP COUNTER* 

(5) MAXIMUM NO* OF VARIABLES TO BE SOLVED* 

NO* OF DIFF. EONS. INITIALLY NOT SOLVED. 

E.G. -2 = DELAY INT* TKE AND DISS* 

UNTIL C4EDSU IS SATISFIED* 

WHEN NCNADD IS SATISFIED, NEQADD » NEQAV2 
»«“!* SATISFIED. NEOAdS . SeQAvI t , 

<1) NO* OF DEP* VAR. TO BE INTEG. IN QKNUIN 
(5) NO* OF SCALARS TO PRINT ACROSS A PAGE (MAX * 8) 

DO NOT USE MIXING LENGTH THEORY FOR DIFF* COEF * 

USE MLT FOR SOLUTION OF DIFF. COT. 

DELAY USING MLT UNTIL E1E2SU IS SATISFIED. 

NO* OF 'NB' BYTE WORDS IN TITLE FOR EACH tifp u&d 

rfilf*r AT WHICH MAXIMUM ITERATIONS TOOK PLACE* * * 

CURRENT PASS (FOR PRINT POINT) * 

STARTING LOC* IN DEP* VAR* MATRIX FOR THIS uartabi r 
<1> IMPLICIT INTEGRATION SCHEME BEING WED. 


NEQAV3* 

NE0AV3* 


0 

1 


I- (4) 


■ * »■» 
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OF POC.i V4UMU1 * 


ALPHANUMERIC I ARRAY ENTRIES 


IARRAM 

ENTRM 


NAHE 


DESCRIPTION 


316 


NIT 


302 

94 

323 

88 

191 


NITER 

NIZS 

NLAST 

NLINE 

NM 


306 

60 

206 


NMBJAC 

NMBOUT 

NMDL 


190 


NMOUT 


193 

16 

340 

55 

233 

19 

325 

142 

40 

38 


NM2 

NNODE 

NNROU 

NODE 

NODES 

NOE 

NOUEDG 

NOUTPR 

NOUTS 

NOUTVC 


- 1 * FLAG FOR DEBUG PRINT OFJ 
A* JACOBIAN 

B. L/U DECOMPOSITION 

- ITERATION COUNTER. 

" <200) NO. OF IZ ENTRY POINTS THAT CAN BE STORED. 

- KOIE FOR LAST STATION UHEN DX .LT. H. 

- (60) LINE COUNT OUTPUT CONTROL. 

- <3> TYPE OF ELEMENTS IN SOLUTION. 

- 2 = LINE (ONE-DIMENSIONAL) ♦ 

- 3 - TRIANGLE (TWO-DIMENSIONAL). 

- (NBANDftNNODE) NO. OF ENTRIES IN JACOBIAN. 

- (30) NO. OF VARIABLES TO BE PRINTED. 

- (8) ALLOW EXTRA STORAGE IN IZ(71) AND IZ<72-76> 

LENGTH OF IZ(71) = MAXIMUM ( NODE t KODEtNMDL) 

IZ<71 ) +2 * NODE THRU IZ<71> + 7 * NODE 
CONTAIN THE REYNOLD STRESSES. 

LENGTH OF IZ< 72-76) = MAXIMUM ( NODE » < NODE/2 )*NMDL> 

- (3) * PRINT OUTPUT IN GEOMETRY FORM. 

2 = PRINT OUTPUT IN NODE NO. SEQUENCE. 

- NM**2. USED FOR STORING FULL MATRICES. 

- NUMBER OF NODES IN SOLUTION. 

- NO. OF COLS. (ROWS) ON ONE SIDE OF JUNCTURE CORNER. 

- (100) VARIABLE DIMENSIONING PARAMETER IN FEDIMN. 

- (NODE) VECTOR LENGTH FOR REFINE GRID GENERATOR. 

- NO. OF EQUATIONS BEING SOLVED FOR DEP. VAR. ' NP ' ♦ 

- FLAG TO STORE UEDGE FROM CP INFORMATION INTO U1 VECTOR 

- (100) NO. OF SCALARS TO PRINT IN OUTPUT. 

- (10) NO. OF OUTPUT VECTORS TO PROCESS AT ONE TIME. 

- (8) NO. OF OUTPUT SCALARS ACROSS PAGE AT PRINT STATION 

WHEN OUTVEC IS CALLED FOR NMOUT = 2. 


♦ 


30 

157 

336 

198 

199 
20 

153 
161 
371 
313 
328 

67 

154 
252 
146 


NP 

NPASS 

NPDBUG 

NPGRDT 

NPGRDV 

NPRNT 

NPUNCH 

NPVSX 

NPVSXT 

NR 

NRJACB 

NS 

NSD 

NSELEM 

NSFDBE 

NSKIP 


- DEP. VARIABLE BEING SOLVED AT THIS TIME. 

-NO. OF PASSES THRU DERVBL. 

- NO. OF ITERATIONS FOR DEBUG IN DERVBL AT DEBUG POINT. 

- (4) STARTUP COUNTER USED IN PRSGRD. 

- (4) STARTUP COUNTER USED IN PRSGRD, 

- (132) NO. OF PRINT POSITIONS ON A LINE OF OUTPUT. 

- SET » 7 IF ELEMENTS AND NODES ARE TO BE PUNCHED IN DIMEN. 

- <2) NO. OF PRESSURES IN P VS X TABLE. 

- NO. OF ENTRIES IN EACH TABLE. (USED WITH NTABPT ) ♦ 

-NO. OF PRINTS FOR IMPLICIT INTEGRATION DEBUG. 

- (1) COMPUTE JACOBIAN EACH ITERATION. 

- GENERAL DUMMY PARAMETER. 

- INITIALIZATION CODE IN DFCFBL. 

- NO. OF SUPER ELEMENTS IN GRID GENERATOR. 

- RESET CONDITION FLAG IN 'FINDBE' » 

- (NODE) NO. OF BOUNDARY LOC. / DEP. VAR. 


47 


ORUW~ 

OF POOR QUALITY 


ALPHANUMERIC I ARRAY ENTRIES 


iarrayT 

ENTRY1 


NAME 


DESCRIPTION 


64 

2S1 

18? 

155 

394 

140 

203 

62 


ASH 

4SN0DE 

*ISTD 

AS2 

NTABPT 

NTCNTS 

NTCRDH 

WTITL 


(10) STOP PROGRAM IF ANY OUTPUT EXP. IS .GT. NSM. 

NO. OF SUPER NODES IN GRID GENERATOR. 

(8) NO. OF STANDARD MATRICES TO BE STORED. 

INITIALIZATION CODE IN DFCFBL. 

NO. OF TABLES TO BE READ IN CPSTUP. 

STARTUP PARAMETER IN CONTES. 

STARTUP PARAMETER IN XYCRDM. 

(10) NO. OF TITLE CARDS TO BE READ IN AND PRINTED AT THE 
BEGINNING OF EACH OUTPUT SET. 


197 WTPRNT 
177 NU2P0S 


178 


NU3P0S 


9999? * DO NOT PRINT INTEGRAL PARAMETERS IN TRBTHK . 
(20) MAX. NO. OF ENTRIES FOR VAR. GEOMETRY DEFINITION. 
(CROSS PLANE) 

(20) MAX. NO. OF ENTRIES FOR VAR. GEOMETRY DEFINITION. 
(TRANSVERSE PLANE) 


260 

90 

91 
192 
166 


NVAR 

NYY 

NZZ 

N2M 

N3DPNS 


NUMBER OF VARIABLES TO BE DISTRIBUTED OVER REFINED GRID. 
(4) SETS (NODE) / DEP. VAR. LOCATION IN IYY VECTOR. 

(4) SETS (NODE) / DEP. VAR. LOCATION IN IZZ VECTOR. 
(NM*2) USED FOR STORING SYMMETRIC MATRICES. 

SET * 1 WHEN JV .LE. (NEQKNN+NEQADD) IN IMPLCT. 


48 


rrr.'T- C3 

q[. pOGR QUALITY 


REAL SCALAR VARIABLES (RARRAY) 
CROSS REFERENCE LIST 


* * * NOTE * * * 

VALUES ZN PARENTHESIS ARE DEFAULT CONDITIONS* 


ALPHANUMERIC RARRAY ENTRIES 


RARRAM 

ENTRlj 


NAME 


DESCRIPTION 


156 

5 

3 

393 

275 

86 

137 

27? 

20 ? 

176 

311 


AINF 

AJ 

ALC 

AM8 

AOMGEX 

AVD 

AVDP 

BEXP 

IBLTH 

jCBTOKJ 

Icci 


- REFERENCE SPEED OF SOUND* 

- (778*28) JOULES CONSTANT. 

- (MIN* SIDE) CHARACTERISTIC ELEMENT SIZE* 

- TERM 8 MULTIPLIER IN PPRES FOR PP EQUATION* 

- (2.0) EXPONENT ON WALL DAMPING FACTOR 'OMEGA' 

- (25.3) DAMPING FACTOR IN DFCFBL* 

AVDP * AVD * SORT ( RHO/RHOW AL ) « ANULOC / ENUT 

- (4.0) GAMMA * 0*01 * DELTAY * ( Y/DELTA>**BEXP 

- BOUNDARY LAYER THICKNESS* DELTA* 

- (4.184) SPECIFIC HEAT BRITISH TO MKS* 

- COEF* FOR REYNOLD STRESSES COMPUTED FROM CPHI1 AND CPH2 


312 

313 

314 
365 
356 

303 
302 
322 

304 
364 
211 

83 


CC2 

CC3 

CC4 

CD 

CE 

CHIEFS 

CHISTP 

CH1TST 

CIMPTH 

CK 

CF0V2 

COMPX 


84 


COMPY 


124 

70 

62 


CON 

CONRHO 

:onv 


IN DIMEN* 

SAME AS CC1. 

SAME AS CC1* 

SAME AS CC1* 

(0.0?) TKE - DISS* COEF. 

(1*0) DIFF. COEF. MULTIPLIER FOR TKE. 

(1.0E-4) CONVERGENCE FACTOR FOR IMPLICIT INTEGRATION* 

(4.0) MAX. NO. OF ITERATIONS FOR INCREASING STEP SIZE. 

(10.0) MAX. NO. OF ITER. BEFORE DECREASING STEP SIZE. 
(0.5) RELAXATION FACTOR FOR IMPLICIT INTEGRATION. 

(1.0) TKE - DISS. COEF. 

SKIN FRICTION 

COMPRESSION FACTOR FOR OUTPUT COL. VECTOR 

INDICATES PERCENT OF X3 AXIS TO BE USED TO SHORTEN 
SPACING INTERVALS. 

COMPRESSION FACTOR FOR OUTPUT ROW VECTOR. 

SAME AS COMPX* BUT FOR X2 AXIS* 

(0*435) KARMANN'S CONSTANT USED IN MLT IN DFCFBL. 

IF »GT« 0.0* SET ALL RHO - CONRHO. 

(1.0) OUTPUT SCALE FACTOR * 1.0 / REFL. 


ORIGINAL PAGE IS 

OF POOR QUALITY 


ALPHANUMERIC R ARRAY ENTRIES 


DESCRIPTION 


REAL NO# EQUIVALENT OF NCONV (NO# OF ITERATIONS) • 
(2.8) EMPIRICAL CONSTANT FOR COMPUTING CC1 - CC4 

!°Ai5L EMF * IRICAL constant FOR COMPUTING CC1 - CCA 
(.24*32.174) SPECIFIC HEAT OF AIR. 

(3.445*32.174) SPECIFIC HEAT OF HYDROGEN. 

SPECIFIC HEAT COMPUTED IN CPINIT. 

REFERENCE SPEDIFIC HEAT. 

(4186.0) SPEC. HEAT CONVERSION FACTOR. 

(1.0) ENTHALPY CONVERSION FACTOR. 

(4.725E-4) PRESSURE CONVERSION FACTOR. 

(16.02) DENSITY CONVERSION FACTOR. 

(1.0) TEMPERATURE CONVERSION FACTOR. 

(0.3048) VELOCITY CONVERSION FACTOR. 

(1.0) SCALE FACTOR FOR DISSIPATION DIFF. COEF 
(1.44) DISS. FCT. PRODUCTION TERM COEFFICIENT! 

(1.0) TKE PRODUCTION TERM COEFFICIENT. 

(1.92) DISS. FCT. DISSIPATION TERM COEFFICIENT. 

(1.0) TKE DISSIPATION TERM COEFFICIENT. 

F0R FUSION COEFFICIENTS. 

(30000*0) TKE - DISS. STARTUP POSITION IN DFCFBL 
(DSTART*SSINIT)**3) IN HFCFBL. 

(DELTSmiELMLT) DEN. AND DISS. TERMS FOR TKE. FPS 
ARE SET TO 0.0 WHEN. DELUSQ , LT. DELCHK. 

(1.0E-5) SCALE FACTOR FOR DELCHK. 

(2.0) PERCENT INTERVAL FOR PRINTOUT. 

DISPLACEMENT THICKNESS. 

FIRST DERVBL PASS FOR LARGEST GRAD*U**2 ON ELEMENT 

(101.0) PERCENT OF TD TO BE USED FOR PLOTTING STATIONS 
(5. 0/9.0) DEGREES RANKING TO DEGREES KELVIN 

(10.0) SCALE FACTOR FOR DCHECK. 

SET .GT. 0.0 FOR CONSTANT PRESSURE FIELD. 

(2.3244) ENTHALPY BRITISH TO MKS. 

(-0.01) ABS(EFMULT) BECOMES LARGEST LEVEL 
OF TKE INITIALIZED. 

(UINF**2) TKE NON-D FACTOR. 

DIMENSIONAL ENERGY FOR U1 VELOCITY. 

NON-DIM. ENERGY FOR U1 VELOCITY. 

DIMENSIONAL I ZING FACTOR FOR ENERGY. 

( XMUINF*RE*UINF*UINF/ ( GtALC** ( NM-1 ) ) 

CONVERGENCE FACTOR FOR INTEGRATION STEP. 

UINFtW/ALC) NON-D FACTOR FOR DISSIPATION. 

(1.0E-n»> SCALE FACTOR ON XMUINF FOR MINIMUM 
LEVEL OF DIFFUSION, 


cnicrm przz vs 

OF POUR QUALITY 


ALPHANUMERIC RARRAY ENTRIES 


RARRAY 

ENTRY 

NAME 

DESCRIPTION 

9 5 

EPTEST 

- * EPMULT - ZERO TEST FOR DISSIPATION USED IN DERVBL 

68 

EP4MD 

- <1.0) MULTIPLIER FOR XMDOT IN PRSORD. 

368 

E8CP 

- (3*0) SCALE FACTOR IN TKE AND DISS* LENGTH. 

171 

EULER 

- (PSTAG/(RH0INF*UINF**2> EULER NUMBER. 

145 

E1E2SW 

- (30000.0) STATION AT WHICH TO SET NE1E2 = 3 - NE1E2 

1 

FACT 

- (ALC) NON-DIMEN. FACTOR. 

80 

FACTH 

- (1.0 / (CP0INF*T0F1NF>> 

59 

FACTMU 

- (RHOINF*UINF*ALC> 

79 

FACTP 

- (1.0 / FACTMU) 

163 

FTTOCM 

- (30*48) FEET TO CENTIMETERS. 

162 

FTTOXN 

- (12*0) FEET TO INCHES. 

164 

FTTOMT 

- (0*3048) FEET TO METERS* 

189 

FI 

- Y-COORDINATE OF FI CURVE. 

329 

F10 

- LAST VALUE OF FI CURVE. 

327 

F2 

- Y-COORDXNATE OF F2 CURVE* 

31 

6 

- (1*0) GRAVITATION CONSTANT. 

277 

GAMEXP 

- (9.0) GAMMA ■ 1.0/ (1*0 + 6AMFAC*(Y/DELTA)**GAMEXP> 

276 

6AMFAC 

- (1.0E-20) SEE 6AMEXP DEFINITION. 

60 

6AMMAF 

- (1*4) FACTOR USED IN GAS LAWS* 

263 

QSCALE 

- (0*1) SCALE FACTOR FOR U1 MULT* OF GRADIENT (PHI) IN 

353 

6UMULT 

- NODE TO ADD GRADIENT (PHI ) TO RHS OF U2r U3 EONS. 

142 

61 

- Z-COORDINATE OF G1 CURVE. 

381 

610 

- LAST VALUE OF G1 . 

328 

62 

- Z-COORDINATE OF G1 CURVE* 

187 

622 

- VARIABLE GEOMETRY FACTOR, 

188 

623 

- VARIABLE GEOMETRY SCALE FACTOR. 

140 

632 

- VARIABLE FACTOR. 

141 

633 

- VARIABLE FACTOR. 

15 

H 

- CURRENT TRIAL STEP SIZE. 

16 

HMAX 

- (2.0) MAX. PERCENT OF TD TO USE AS STEP SIZE. 

45 

HS 

- CURRENT STEP SIZE. 

7 

HSXNXT 

- (1.0E-5) START INTEGRATION STEP SIZE AT THIS VALUE. 

12 

HT 

- OUTPUT VAR. FOR TIME STEP ■ HS * FACT / REFL 

186 

H21 

- (1.0) GRID GROWTH SCALE FACTOR. 

386 

H21L 

- LAST VALUE OF H21. ( USED IN XYCRDM ) 

139 

H31 

- (1.0) GRID GROUTH SCALE FACTOR. 

387 

H31L 

- LAST VALUE OF H31. ( USED IN XYCRDM > 

272 

0ME6XP 

- (1.5) TKEEXP ■ 2.0 * < 2.0 - OMEGXP > 

2 

ONE 

- (1.0) PROGRAM CONSTANT. 

396 

080 

- ADD STRESSES TO PP EQUATION IN PPRES. 

375 

08H18Q 

- (1.0 / H21**2> 


PPRES . 
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ALPHANUMERIC R ARRAY ENTRIES 


RARRAY 

ENTRY 


NAME 


DESCRIPTION 


136 

OSMAX 


271 

OSUSQ 


362 

0S12 

- 

361 

0S6 

■ 

363 

0S60 

- 

339 

PCFACT 


174 

PDFTOC 

- 

170 

PDFTOK 

- 

36 

PEDD1M 

- 

39 

PEDGE 

-■ 

268 

PHICOD 

- 

82 

PIBAR 

- 

9 

PINF 

- 

180 

PMSKGS 

- 

340 

PPFACT 

- 

99 

PPRCON 

- 

100 

PPRIME 

- 

67 

PR 

- 

185 

PRDIS 

- 

166 

PSFTOA 

- 

169 

PSFTOI 

- 

160 

PSFTON 

- 

167 

PSFTOT 

- 

351 

PSTAG 

- 

20 

PTII1 

- 

345 

PUMULT 

- 

179 

RADCON 

- 

21 

RE 

- 

43 

REFL 

- 

47 

REFLRE 

- 

285 

RELAX 

• 

289 

RHOIM 

• 

10 

RHOINF 

- 

157 

RHOUIN 

• 

360 

RHSCAL - 

» 

199 

RNULOC - 

» 

290 1 

ROMULT ■ 


119 1 

ROUALC * 


116 1 

RR 


32 

v’TCONl ■ 


56 

RTCONS - 

1 

117 

iTOHMl - 

1 


■ ADD PC TERM TO U2t U3 EQUATIONS. 

• <0.01602) POUNDS/FT *$3 TO GRAMS/CM**3 

• <16.02) P0UNDS/FT**3 TO KG/M**3 
DIMENSIONAL PRESSURE * PEDGE * PST AG. 

NON-DIM. PRESSURE AT PRESENT STATION, 

U.O/H) GRADIENT PHI MULTIPLIER AT END OF PPRES. 

UTAU COEF. FOR COMP. U1 IN JNCXNP. 

FREESTREAM PRESSURE. DEF, = 1ST VALUE IN P VS X TARt F 
<1.0/2, 2) POUNDS / KG. TABLE 

USE PP TERM IN U2» U3 EQUATIONS. 

<RH0INF*UINF**2/ALC> 

PRSSURE GRADIENT COMPUTED IN PRESSURE ROUTINE. 

(1.0) PRANDTL NUMBER. 

<1»3) DISSIPATION PRANDTL NUMBER. 

<4.725E-4> P0UNDS/FT**2 TO PSIA 
<6.924E-3> P0UNDS/FT442 TO P0UNDS/IN**2 
<47,88) P0UNDS/FT**2 TO NEWT0NS/M*»2 
<0.3591) P0UNDS/FT**2 TO TORR, 

< PINF+O. 5#RH0INF*UINF**2 > 

PRINT TIME PARAMETER IN 'IMPLCT'. 

WHEN = 1» USE EDDY VISCOSITY FOR U2r U3 DIFF. TERM 
<57.2957759) CONVERSION FACTOR RADIANS TO DEGREES 
<RHOINF*UINF*ALC/XMUINF> REYNOLD'S NO. EES * 

<1.0) REFERENCE LENGTH. 

<RH0INF*UINF4REFL/XMUXNF) REYNOLD'S NO. BASED ON REFL 
0.4) RELAXATION FACTOR FOR GRIB GROWTH IN LOOKAW 

(Finf,xha/?runiu*tatnc 4 TKE ’ EPS “ons! N in°Serwl. 
rhoinf X * A uinf MIW * TSINF> freestream oensity. 

??S” - U8| R r0CAL R V ?M2si?Y 8 ?SR S VAN DRIEST^BANFING FACTOR 

RH«N?T^Nrr^c«i u * ToFiNF > »s: 

CPH / CPA 
2.0 * G * AJ 

UINF**2 / <RTC0N1*CP0INF*T0FINF > 

RR * < TOH/TOA - 1.0 ) 


onicsria f/»c 4 s is 

of poor quality 

ALPHANUMERIC RARRAY ENTRIES 


4 


R ARRAY 


NAME 


DESCRIPTION 


ENTR^ 

280 


RUEDSW 


(10,0) IF < Y/DELTA .GT. RUEDSW > 

EPSDIM = 6AMREF 4UEDGE4DST ARtGAKMA 


138 

28 

129 

210 

301 

190 

50 

75 

73 

74 
72 
71 

321 


RUESQ 
RUN IV 
SCT 

SHPFAC 

SIMPLT 

SLOPE 

SSINIT 

STLCON 

STLDCR 

STLDEX 

STLDTR 

STLDVR 

STPMLT 


- RH0ED6 $ UEDG * UEDG 

-(1545. 33*32. 174) UNIVERSAL GAS CONSTANT. 

- (1.0) CONSTANT SCHMIDT NUMBER. 

- SHAPE FACTOR. 

- (TO) STATION AT WHICH TO START IMPLICIT INT. 

- SLOPE OF VARIABLES COMPUTED IN LOOK. 

- (HSINIT / FACT) 

- WHEN .GT. 0.0» HAVE CONSTANT VISCOSITY. 

- (204.0) REF. CON. TEMP. IN SUTHERLAND. 

- (1.5) EXPONENT USED IN SUTHERLAND* 

- (492.0) REF. TEMP. USED IN SUTHERLAND. 

- (1.163E-5) VISCOSITY USED IN SUTHERLAND. 

- (1.1) STEP SIZE MULTIPLIER WHEN 

ALL VARIABLES CONVERGE ON 1ST ITERATION. 


374 

372 

35 

22 

206 

4 

23 
48 
42 
90 

300 

305 

24 

146 
58 

147 
40 

284 

347 

155 

111 

26 

394 

398 

397 

371 


TAREA 

TBAR 

TD 

TF 

THETA 

THK 

TIME 

TIMESV 

TKED6E 

TKEINF 

TMNTS 

TMULT 

TO 

TOA 

TOFINF 

TOH 

TRATIO 

TSADD 

TSCALE 

TSINF 

TWELVE 

TWOPI 

T2FIX 

T2PFIX 

T3FIX 

UBAR 


TOTAL COMPUTATIONAL AREA. 

MASS WEIGHTED AVERAGE TEMPERATURE. 

(1.0) TOTAL SOLUTION TIME ( DISTANCE ) FROM TO. 
(TF*TO+TD) FINAL STATION. 

MOMENTUM THICKNESS. 

(1.0) DEF. NON-DIM. THICKNESS OF ELEMENTS. 

CURRENT STATION. 

SAVED TIME LOCATION FOR IMPLICIT INTEGRATION. 

LIMIT ON FREESTREAM EDDY VISCOSITY COMP. IN DFCFBL. 
(UINFII2) 

TIME (MINUTES) OF CPU USED. 

(1.04) STEP SIZE MULTIPLIER. 

STARTING TIME (DISTANCE). 

(533.0) AIR REFERENCE TEMP. FOR COMPUTATIONS IN DIMEN. 

(533.0) REFERENCE TEMPERATURE. 

(520.0) H2 REF. TEMPERATURE FOR COMPUTATIONS IN DIMEN. 
1.0 ♦ (GAMMAF-1.0) * XMACHS**2 / 2.0 

INITIAL LEVEL FOR SCALE FACTOR 'TSCALE' 

SCALE FACTOR FOR T2FIX» T3FIX AND T2PFIX IN PPRES. 
STATIC TEMPERATURE COMPUTED IN CPINIT. 

(12.0) LENGTH SCALE USED IN BRDSHW. 

(2.0BPI) 

ADD U2 CONVECTION TO RHS OF PP EQUATION. 

ADD U2' AND U3' TO RHS OF PP EQUATION. 

ADD U3 CONVECTION TO RHS OF PP EQUATION. 

MASS WEIGHTED AVERAGE VELOCITY. 
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RARRAlf NAME 


ENTR 


385 

UCMULT 

203 

UED 

63 

UEDOE 

27 

UINF 

354 

UINFX 

202 

MUALL 

458 

U1MIN 

282 

U2STRS 

120 

VARB 

348 

VCMULT 

104 

YELCST 

177 

JLBTON 

178 

/LBTOP 

346 

SLDMLT 

102 l 

JSTART 

383 L 

ISMAX 

273 S 

:delta ■ 


EPS, 


125 

109 
61 

154 

373 

66 

172 

110 
38 
98 

52 
11 

201 

330 

324 
198 

53 
382 

325 


LAM 
MA 

MACHO 
MACHS 
MDOTC 
MF 

MFACT 

MH 

MU INF 
PRIME 
SCALE 
T 

TC 

MULT 

NRMAD 

PLUS 

SCALE 

MULT 

NRMAD 


.ALPHANUMERIC RARRAV FNTfrrrc 

description 

' w’wSSffJ? 0 " TERM U». TKEi 

- FREESTREA« F UtLOCI T r N ® RB8HU ’ 

- (UINF) 

■ VELOCITY just of wall. 

■ AD^XoL^STRESSES^O^f U^EOuIf In'* ' IS CALlEI ' 

s NS - 

E ™ s - 

sas sS» «• 

“ ?eS 0 St\ D Ht T S„^h g §?“?I I0NS - 
«»th a H 

‘••“’Htfraia -ggS’Jtr 1 ""- 

<0.09) CONSTANT USED JW^FCFBL^ 

MA 2 C 8 H 9 N 7 iMB M E°R L ? CULAR WEI6HT » *'«< 

J-fJpAL MACH NUMBER. 

AVERAGE MAS FLOW. 

(UlN?2c22i ECUUAR weight OF FLUID. 

SsSpwfap” 

aSSSSSssr - 

SaS*:r»'- 
s.raa sss swsF'~ wa: 

(010) MSPLACENeSt OF*mmS MULT »f>-IER. 

1SFLACEHENT OF COORD. IN Z DIRECTION. 
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THE FOLLOWING IS A LIST OF ENTRIES AT THE BEGINNING OF THE 
12 ARRAY (1ST 'NIZS' ENTRIES) WHICH CONTAIN ENTRY POINTS IN 
THE REMAINDER OF THE IZ ARRAY (STARTING AT IZ(NIZS41)> 

WHICH ARE THE STARTING LOCATIONS FOR THE VARIABLE LENGTH 
VECTORS USED IN THE PROGRAM. 


* * * NOTE # $ 4 

VALUES IN PARENTHESIS DENOTE DIMENSION OF VECTORS. 


LOCATION IN I ARRAY FOR DIMENSIONALIZING SCALARS. 


113 

170 

194 

191 

206 

371 

62 

112 


- KPLVAR » 


5 *> - 


NBC> 25 
NDP» 31 
NM» 192 
NMDLt 55 
NPVSXT > 27 

NTITL» 176 
NWALLSt 90 


KROWf 50 
NCALLSt 210 
NEQr 58 
N2M» 193 
NODE* 142 
NSKIPf 121 
NTKSt 177 
NYY * 91 


LCOL f 47 
NCOMPG » 59 

NEQKNNt 46 
NM2» 60 
NOUTPR f 162 
NSPECf 189 
NU2F‘0S * 178 
NZZ » 252 


LG» 
NCFTABf 
NF, 
NMBOUT » 
NPVSX» 
NSTD» 
NU3P0S » 
NSELEM 


DEFINITIONS FOR OTHER DIMENSIONS BASED ON ABOVE SCALARS. 

MAXDEP = MAXO ( NODE» N0DE4NMDL ) 

MAXVAR “ MAXO ( NODE » MAXDEP/2 ) 

MGLC = MAXO ( LCOLf LG*NWALLS ) 

NEMD = NODE * 1.81 
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1Z ARRAY ENTRY POINTS 


IZ 

NTRY 


1 

2 

3 

4 

5 


NAME 


I COL 
I ROW 
IFMTHD 
ITITLE 
IIPXNT 


6 j IKBNO 

10 | IINCOL 

11 ilINROW 
14 ' IPSIBD 
18 jlCPTAB 

19 ITTAB 

20 1 I IUSED 


25 


IIBND 


26 1 1 INODE 
27 jlJBOND 
28 ! IKEYCL 

29 IKEYDG 

30 | IKEYRW 
32 | I NUN 
33 I INDEX 
34IIINDRU 

INOCOL 
1 1 ELS 
HELEN 
IIBORD 
IRHOE 


DIM* 


DEFINITION 


(LCOL ) - DISCRETIZATION COLUMN LOCATIONS 
|( KRQW ) I- DISCRETIZATION ROW LOCATIONS. 

(NFtNMB0UT+20) - HEADINGS FOR OUTPUT VARIABLES* 
(NTITL420+20) - TITLE FOR START OF EACH OUTPUT PHASE* 
l<NDP*10+NSPEC> - 1ST NDP LOCATIONS - DEPENDENT VARIABLES* 

- 2ND NDP LOCATIONS - PREDICTED VARIABLES. V "* IAHLES * 

- 3RD NDP LOCATIONS - LOCATION OF IND* VARIABLES. 

- 4TH NDP LOCATIONS - LOCATION OF PREDICTED VARIABLES* 


(NDP ) 
(LC0L*2> 
(KRQW ) 
(NSELEM) 
(NCPTAB) 
(NCPTAB) 
(NODE ) 


35 

36 

37 

38 

39 


43 

45 

46 

47 

48 

49 

50 
52 


IC200 

ISMSTR 

IBCNST 

IDIF 

IYY 

IZZ 

IX1P2 

lYDIn 


NO. OF BOUNDARY NODES / DEP* VARIABLE . 

NO. OF NODES PER COLUMN. 

NO* OF NODES PER ROW. 

* EVERSAL KE * r 0R TYPE II ELEMENTS. 
SPECIFIC HEAT TABLE ENTRIES. 

TEMPERATURE TABLE ENTRIES. 

COUNTER USED IN OUTPUT ROUTINE. 

(NODEtNEQ) - RE-ORDERED NODES / DEP. VAR. TO ACCOUNT 
J - FOR BOUNDARY CONDITIONS. 

!f!!?5* NM 1 ' ARRAY 0F ele « e NT CONNECTIONS (NM/ELEMEN1 ) • 

' U NODE SOLUTION ORDER USED IN BANCHO. 

■ COLUMN KEYS FOR BANCHO. 

- DIAGONAL KEYS FOR BANCHO* 

• ROW KEYS FOR BANCHO. 

• TEMPORARY STORAGE. 

ORDER OF NODES BY COLUMNS FROM LEFT TO RIGHT. 
ORDER OF NODES BY ROWS FROM TOP TO BOTTOM. 

OUTPUT COLUMN POSITION OF NODES BY ROWS. 

NO. OF ELEMENTS CONNECTED TO NODES. 

LIST OF ELEMENTS CONNECTED TO NODES. 

° F B0Rt ' ER NOr ' ES 1N COUNTER-CLOCKWISE ORDER. 
TEMP. STORAGE FOR RHOtKt RH04EPS. 

2*NM + N2M + 2*NM2 + 2*NM**3 + 50 


{(NODE 
(NODE 
(NODE 
(NODE 
(NODE 
(NODE 
(NODE 
(NODE 
(N0DEK2) 
(N0DE*6> 
(NODE > 
|(N0DE*2) 
NSTDMX 


(NSTDMX) j- STORAGE FOR STANDARD MATRICES. 

|( NSTD+5) I- ENTRY POINTS IN IZ FOR STANDARD MATRICES. 

NBCND * NSKIP + 6*NBC + 1 
( NEQ4NBCND+ 10) - STORAGE FOR BOUNDARY CONDITIONS* 
NDIFST = NODE ¥ (NEQ+NEQKNN) + NEQ + 10 

DI ! F UfION COEFFICIENTS / DEP. VARIABLE. 
(NYYtNODE) - 4 SETS OF VALUES / DEP. VARIABLE. 
NZZ*NODE> - 4 SETS OF VALUES / DEP. VARIABLE. 

KNEMD > - ELEMENT LENGTHS COMPUTED IN GEOMFL . 

|( NODE > I- NOT USED. 


56 


IZ ARRAY ENTRY POINTS 


IZ I 

ENTRY INANE 


53 IPCOL 

54 i I PROW 
54 IAMXLT 

41 -ITBRTK 

43 XU2POS 

44 IU2VAL 

45 IU3P0S 
44 IU3VAL 
71 ‘iOUTl 


72 I0UT2 

^ *. 

TO ' . 

74 ;I0UT4 

77 I AREA 

78 .ICP 
7? ;ih 
80 iXPSI 
82 IQ 

IQP 
IRHO 
ITENP 
IRHSP 
ITK 

XXI COR 
IX2C0R 
I PRESS 
XANU 
XNORNY 
INORNZ 
NJST 
IPRVAL 
IOL 
XQPL 
I MEL 
1U 

I PR OR D 
IYNOD 


DIM, 


DEFINITION 


(LCOL ) - COLUMN COORDINATES FOR OUTPUT PAGE, 

(KROW ) - ROW COORDINATES FOR OUTPUT PAGE, 

(N0DE*4) - MIXING LENGTH? DUDY? WALL DAMPING AND DIFFUSION 
LENGTH, 

(MGLC*NTKS+2) |- STORAGE FOR INTEGRAL PARAMETERS. 


83 

84 

85 
84 
88 
8 ? 

90 

91 

92 
95 
94 
99 

100 

101 

102 

103 

104 

105 
108 


109 

110 


(NU2P0S*2> 
(NU2P0S42) 
<NU3P0S*2> 
(NU3P0S*2) 
(MAXDEP) 


(MAXMAR) 


- DOWNSTREAM PLOS, FOR CROSS COORD. CHANGE. 

- SCALE FACTOR FOR CROSS COORD. CHANGE. 

- DOWNSTREAM POS. FOR TRANSVERSE COORD. CHANGE. 

- SCALE FACTOR FOR TRANSVERSE COORDINATE CHANGE. 
TEMPORARY STORAGE. 

<I0UT1+2*N0DE> THRU ( I0UT1+7*N0DE) CONTAIN 
REYNOLD STRESSES. 

TEMPORARY STORAGE. 

( I0UT2 + NODE) CONTAINS GAMMA IN SUB. DFCFBL. 


(MAX VAR) 


(NEMD 
(NODE 
(NODE 
j ( NODE 
I (NODE 
! ( NODE 


(N0DE*3) 


IGE0M1 

IGE0M2 


(NODE 
(NODE 
(NEMD 
(NODE 
(NODE 
(N0DE*2) 
(NODE ) 
( N0DE42) 
(N0DE42) 
(N0DE*3) 


TEMPORARY STORAGE. 

AREA OF ELEMENTS COMPUTED IN GEOMFL. 

NODAL VALUES OF SPECIFIC HEAT. 

NODAL VALUES OF ENTHALPY. 

TEMPORARY STORAGE FOR EQUATION SOLVING IN STRF. 
TEMPORARY STORAGE. 

TEMPORARY STORAGE. 

NODAL VALUES OF DENSITY? RH04U1 AND AVE. RHQ*U1. 
NODAL VALUES OF TEMPERATURE. 

RIGHT HAND SIDE FOR EQUATION SOLVING IN STRF. 
ELEMENT THICKNESS DISTRIBUTION. 

NODAL VALUES OF TRANSVERSE COORDINATES. 

NODAL VALUES OF NORMAL COORDINATES. 

NODAL VALUES OF PRESSURE. 

NODAL VALUES OF LAMINAR VISCOSITY. 

NORMALIZED TRANSVERSE COORDINATES. 

NORMALIZED VALUE OF CROSS COORDINATES. 

STORAGE FOR RHS IN STRF. 


(NCOMPGt 10) |- RESTART STORAGE FOR PRSGRD VARIABLES. 
(N0DE48) - STORAGE FOR DERVDX ROUTINE. 

(NODE*?) - STORAGE FOR DERVDX ROUTINE. 

(NODE ) - TEMPORARY STORAGE FOR RHO * U2. 

(NODE ) - TEMPORARY STORAGE FOR RHO * U3. 

(N0DE*3) - NODAL VALUES OF DPDX. 

(NODEtNWALLS) |- TRANSVERSE COORD. USED IN CONTES? DFCFBL? 
TRBTHK? WLFLX8? ETC. 

(NEMDtNM) I- NATURAL COORDINATE DERIVATIVE COMP. IN GEOMFL. 
(NEMDtNM) - NATURAL COORDINATE DERIVATIVE COMP. IN GEOMFL. 


XZ ARRAY ENTRY POINTS 
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IZ 

k 

:ntry 

NAME 

hi 

I PLOTS 

112 

IPLTYP 

113 

IPLSCL 

114 

ISCHMT 

121 

ICALL 

123 

IOMULT 

124 

IOSAVE 

125 

INOUT 

127 

IICNCL 

128 

IICNDX 

131 

IIPAR 

132 

10PAR 

133 

ISUTLD 

134 

IPR 

135 

IMPAR 

136 

I EPS 

137 

I8KNFR 

138 

I8TN 

139 

IX3ST 

140 

IPVSX 

143 

NVDEP 

146 

IMXLT 

150 

I8SAV 

151 

IPSAV 

197 

IDOUTP 


DIM. 


DEFINITION 


(KPLVAR*2>|- LIST OF VARIABLES TO BE PLOTTED. 


(KPLVAR*2) TYPE OF PLOT TO BE GENERATED. 

(KPLVARtlO) I- PLOT SCALE FACTORS. 

(NODE ) f- NODAL VALUES OF SCHMIDT NUMBERS. 

(NCALLSX2) - LIST OF LINK NOS. AND ENTRIES TO CALL AT 
I END OF IMPLCT . 

NMB * NMBOUT + 4 

<NMB*2 ) (- LIST OF MULTIPLIERS FOR OUTPUT VARIABLES. 

(NMB ) |- LIST Or VARIABLES TO BE PRINTED IN OUTPUT. 

MM * MAX0(15.NSP£C»NPVSXT+2) 

MMNODE = NODE * MI NO ( MM t NMBOUT) 

( MMNODE ) (- TEMP. STORAGE FOR OUTPUT VAR. AND SOURCE DATA. 
(MGLCB2) r NO. OF NODES / COL. USED IN CONTES. DFCFBL. 

TRBTHK. ETC. 

( NODEBNW ALLS+LCOL ) |- LIST OF NODES / COL. USED IN CONTES. 

DFCFBL. TRBTHK . ETC. 

(NOUTPR) I- LIST OF PARAMETERS TO PRINT AT START OF OUTPUT. 
(N0UTPR44+20) I- TITLE INFORMATION FOR PARAMETERS AT 
BEGINNING OF OUTPUT. 

(10 ) |- LIST OF CONSTANTS USED IN SUTHLD. 

(NODE ) - NODAL VALUES OF PRANDTL NUMBER. 

(NOUTPR) - LIST OF MULT. FOR PARAM* AT START OF OUTPUT. 
(NODE ) I- NODAL VALUES OF TURB. VISC. COMPUTED IN DFCFBL. 
(MGLCB4+10) - SKIN FRICTION DIST. BY COLUMN. 

(MGLC ) [- STANTON NO. DIST. BY COLUMN. 

(NPVSX ) I- LIST OF DOWNSTREAM STATIONS IN PRESSURE TABLE. 
NPNONP * NPVSX + IPCF1TBN0DE 

(NPNONP) I- LIST OF DOWNSTREAM PRESSURES IN F'RESSURE TABLE. 
NDNENM = NODE 4 NEQ*(4*NM+20) 

( NDNENM ) h ELEMENT ADDRESS LOCATIONS USED IN BCONDT. DERVBL 
AND PPRES ♦ 

(LG* 4 > h MIXING LENGTH FROM GETFSL* 

(NEMDtNM) f- TERM 8 STORAGE IN PPRES. 

(N0DEB3) h PHI. U2RHS. U3RHS TEMP. STORAGE IN PPRES. 

(REST OF STORAGE) I - TEMPORARY MATRIX STORAGE. 


V 

( r n • r • < 

■' • • 4 - L*> 

Oh l'o« QUALITY 

IZ ARRAY ENTRY POINTS 


t 


201-240 ARE TEMPORARY LOCATIONS FOR TYPE II GRID GENERATION. 

'NIZS' IN NAME01 NEEDS TO BE SET TO 250 TO ALLOCATE STORAGE 


♦ 


LOCATION IN I ARRAY FOR SCALARS. ] 

261 ~ IBC» 

27 - NSKIP, 55 - NODE t 252 - NSELEM » 

25? - NSIDE» 

251 - NSNODEf 260 - NVAR 

DEFINITIONS 

FOR 

OTHER DIMESIONS BASED ON ABOVE SCALARS i 

KSEL 

ss 

NSELEM + 5 

KSLNSD 

* 

KSEL * NSIDE * 4 

KSNODE 

s 

NSNODE + 5 

KSNVAR 

« 

KSNODE * NVAR 

NNVAR 

» 

(NODE+5) * NVAR 

NODVAR 

s 

NODE * NVAR 


IZ 

ENTRY 

NAME 

DIM. 

DEFINITION 

201 

NZONE 

(KSEL ) 

- SUPER ELEMENT NUMBERS. > 

202 

I SHAPE 

(KSEL ) 

- ORDER OF GENERATED ELEMENTS IN REFINE* | 

203 

IDLZ 

(KSEL > 

- TYPE OF ELEMENTS DESIRED / SUPER ELEMENT. 

206 

NZ 

(KSEL ) 

- NO. OF DIVISIONS IN 3RD DIRECTION. 1 

207 

ITYPE 

( KSEL > 

- TYPE OF SUPER ELEMENT 1 = TRIANGLE. , 




i 2 = QUADRILATERAL. 

208 

NCOL 

(KSEL > 

- NUMBER OF COLUMNS IN SUPER ELEMENT. ' 

20? 

MELEM 

(KSEL ) 

- SUPER ELEMENT TO ELEMENT. 

210 

NELK 

(KSEL ) 

- TYPE OF SUPER ELEMENT LINK. 

211 

G 

(KSEL ) 

- SUPER NODE TO ELEMENT. 

212 

KQ 

(KSEL ) 

- G KEY. 

213 

KF 

(KSEL ) 

- F KEY. 

214 

KTK 

(KSEL ) 

- THK KEY. 

215 

IBCTYP 

(50*IBC) 

- BOUNDARY TYPE KEY. > 

216 

BOUND i 

1504IBC) 

- BOUNDARY COEFFICIENTS. J 

217 

IDUM 

(NODE+5) 

- DUMMY NODE COUNTER IN REFINE. 



NSKIBC = 

* NSKIPtIBC j 

218 

BC 

(NSKIBC) 

- GENERATED BOUNDARY CONDITION ARRAY. 

21? 

NODF 

(NODE > 

- VECTOR OF BOUNDARY NODES. ! 

220 

KTFIX 

(NODE > 

- FIXED DEP. VAR. NODE KEYS. 1 

221 

IRNSTB 

(KSEL48) 

- NOT USED 

222 

JELF 

(NODE ) 

- VECTOR OF BOUNDARY ELEMENTS* 

223 

KBND 

(NODE > 

- SUPER ELEMENT BOUNDARY KEYS. ! 

224 

IRTH I 

(NODE ) 

- NOT USED . 


< 




. 


59 




IZ ARRAY ENTRY POINTS 


ORIGINAL IV.C£ U 

OF POOR QUALITY 


IZ 

ENTRY 

* 

NAME 

DIM. 

22S 

AA 

(KSNODE) 

226 

BB 

< KSNODE > 

227 

SNODE 

(KSNODE) 

228 

SGRID 

(KSNVAR) 

22? 

DUM 

< NNVAR ) 

230 

ELNK 

(KSLNSD) 

231 

MATL 

(NEMD > 

232 

COORD 

(NODVAR) 

233 

INX 

(KSEL > 

234 

INY 

(KSEL ) 

235 

NDVEC 

(NODE > 

236 


(KSNVAR) 

237 

ELCORD 

(NNVAR ) 

238 

KYND 

(NNVAR ) 

23? 

I DM2 

- TEMPORA 

240 

IELNOD 

(NODE )| 


2 . 


DEFINITION 

- non-rect«ngular CARTESIAN COORDINATES mb 1 
NON-RECT ANGULAR CARTESIAN COORDINATES MR* 

- SUPER NODE NUMBER VECTOR. 8 tI8 ' 

VARIABLE ARRAYS RETURNED FROM REETNr 

- DUMMY STORAGE IN REFINE . REFINE* 

SUPER NODE CONNECTION TABLE* 

- Fimre T i?p2r«? R ^f MENT T0 ELEMENT 

ELEMENT node coordinates* 

- NO* OF DIVISIONS IN 1ST DIRECTION. 

2.°cs4^ DIVISIONS IN 2nd direction* 

MnT M M2Jr. C ° NNECTI0NS F0R 6RID GENERATION. 

~ NUT USED 

" ESS5 D ‘ F °* ? R ? D generation* 











ORIGINS. MCE n 
OF POOR QUALITY 



Figure 2. BDINPT flow chart 










Figure 2. BDINPT flow chart (cont.) 




















/ 


^ CALORD ) 


^ COMOC ^ 



EXPNT =0.5 

IF (VSA.VE.EO. POWER )EXPNT*POWER 


NM0D:0 


DO FOR I = 1, N 
RZ(IZ1)=RZ(IZ2)**EXPN1 




RARRAY ( N1 )=RARRAY ( N2 ) **EXPiNT 

. A 

RARRAY(NBRS-l) = 


RARRAY ( NBRS ) **EXPNT 

r-* 


[^40g] 


Figure 2. 3D1NPT flow chart (cont.) 








Figure 2. BDINPT flow chart (cont.) 
































BDIWT: 10 of 10 



Figure 2. BOINPT flow chart (concl.) 









F POOR QUALITY 



Figure 3. OERVBL flow chart (cont.) 








w VI 


OF POOR QUAU1 ‘ 


l30^j 

Compute scalar averages 

§ 

Store diffusion term multipliers. 

Zero out element martices and scalars. 


DO 500 ML * 

1, NM 1 

i 

t 

o out element vector entries 

t 

* 

DO 400 MK 

*— — — —— H 

« 1, NM ] 


Compute $6200? ’1QJ for necessary 

variables. 

Compute |Q$B3OOO}fQ0for gen. & di : 
terms ofTKE and EPS. 


Compute sourc e te rm for U1 mom, eqr 


Assemble diffusion & convection 
. element ve ctors. 


Figure 3. DERVBL flow chart (cent.) 





***** • » uu ♦ 


T U I 


omc-rw- pwz u 

OF POOR QUALITY 


Compute DELLISQ for gen. terms of TKE & EPS 


Compute ENERGY for dependent var. 


Compute elem. contr. to U\ T. M 


Check on magnitude of DELUSQ. 

On first pass only - - compute the 
largest value of DELUSQ and store 
in DELTST. (on completion of 1st 
pass set DELCHK = DELTST * DELMLT ) 


DO 525 ML « 1, NM 


Assemble C Q, Diffusion and 
Convection matrices for Jacobian 


L_D0 450 NP = 1, NEQKNN 

— — £ 

Assemble common te rms in Q RHS 


Figure 3. DERVBL flow chart (cont.) 










DFCFBL: 1 of 6 


Ut ('uun Qw-tl.il Y 



Figure 4. DFCFBL flow chart 
















DO 300 LL«1 , M 


UFCFoL: 3 of 6 


Extract node nos. and nodal var. 
as needed for computations. 
Compute Y+, U+, 1, t, Y/6 , u 
Compute Dim. Eddy Vise. - EPSD1M 



Limit EPSDIM to 
EPSMAX * XDELTA 


Multiply EPSDIM by interim’ ttancy 
factor y . Compute UEDSTR=GAMREF * Ue * 6* *y 
If (Y/ .gt.RUEDSW) set EPSDIM = UEDSTR 


NE1E2: 1 


Y+iYPLUS 


GAMMA 

TKEDGE 


ITKE: 1 


YPLUS: 


Figure 4. DFCFBL flow chart (cont.) 









■- wx*i 13 

OF POOR QUALITY 


uruuL: 4 of 6 






ORJQ , ?J/ l t pr^z (3 
OF FOOh QUALITY 



Figure 4. OFCFBL flow chart (cont.) 


















OP QUALITY 



Figure 5. IMPSLV flow chart 
























/ir , ?o? rr r 1 

or t uuu 


c. ' 
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r, r » 

t -4 


vJAiJ'l Y 


F KQGFv'AM JNCTIJR( JNPU7 f OUTPUT f TAPE 5 a INPUT f TAF'E6*0UTF‘UT ) 


---c-o-m-o-c 

DIMENSION RZ(1)f L ( 400 > 

EQUIVALENCE ( RZ(1)f IZ(1)f 1.(1) ) 

COMMON / VARBLE / I ARRAY ( 00500 ) f RARRAY(OSOO) 

EQUIVALENCE ( I ARRAY ( 00092) f IZSIZF. ) 

COMMON / ARRAYS / IZ( 350000) 

NZSIZE » 350000 

CALL RESET I ( 500 f I ARRAY f 0 ) 

CALL RESET ( 500 f RARRAY f 0,0 ) 

IZSIZE <= NZSIZE 

CALL RESET I ( IZSIZEf IZf 0 ) 

WRITE ( 6t 9600 ) IZSIZE 

9600 FORMAT < 1H f 35H WING JUNCTURE TEST CASE* IZSIZE = , 17 ) 
CALL BDINPT 
STOP 
END 

BLOCK BATA “ 

DIMENSION TITLE( /; f VNAME(3f10)f I NAME ( 6) 

COMMON / F’LOTFE / NTITLf I NAME f NOBDNOf NSf NVf VNAMEf TITLE 
COMMON / CONF’PR / DERV2<10> 

COMMON / AFYSAV / DFYSAVU02) 

COMMON / FSHAPE / D(60) 

COMMON / BISCRT / IDSC(IO) 

COMMON / MANOUT / MANY (20) 

COMMON / NEWGOM / DGEOM(IO) 

COMMON / NPINTQ / JJ(30) 

COMMON / JADRES / JAB (30) 

COMMON / OUTUUM / 

1 FMTl(lO), FMT2(09 ) t FMT3(17)f RNMB ( 1 1 ) » DUM(28> 

DATA FMT1 /4H(1H f4Hf f4H 8f 4H A1 f 4Hf4Xff 4H 1, 

* 4H 2. 4H Ot 4H AIf 4H ) / 

DATA FMT2 / 4H(1H0f 4Hf1Xff 4H1HEff 4H I2ff 4H If 4H 2f 

* 4H 8* 4H AIf 4H ) / 

DATA FMT3 / 4H( 1f4H 4f4HXf f4H »4H lr4H( t4H 4 
1 AH Af4H 4 f 4H f 1HE f 4Hf I3f f 4H f4H 1f4HX ))f3*1H / 

DATA RNMB / 4H Of 4H If 4H 2f 4H 3f 4H 4f 4H 5f 

% 4H 6 f 4H 7 f 4H 8f 4H 9f 4H 10 / 

DATA DUM / 28*0* ) / 

APPENDIX A. 

SUBROUTINES FOR 3DPNS 
TEST CASE 


C" -t> w 


OUlGvtSi'A. IS 

OF POOR QUALITY 


DATA NTITLf I NAME f NOBDNQf N3» NVf TITLEf VNAME / 

1 6, 76f 21 t 27 1 38 t -18?f Gt 0» Of Of 

2 4H ST f 4HANDAf 4HRD P* 4HPINT f 4H VAR t 4H1ABL t 4HES f 

3*1H f 4H f 4HSTAT i 4H1QN f 4HREYNf 4HOLDSf 4H NO.f 

4H RF.f 4HF . Vf 4HEI. . f 4H REF f 4H. VI t 4HSC. f 4MVARI f 

4HABLEf 4H NO.f 12*1H / 

DATA MANY / 39f 47f 56f 91 f 66t 101 f 102f 105f 12*0 / 

DATA DFYSAV / 102*0.0 / 

Df TA DERV2 / 10*0.0 / 

DATA D / 5*0. Ot 0.25f O.Of l.Of O.Of 0.5f O.Of -l.Of 0.075t 
1 9*0. Of 1 .2 f 37*0.0 / 

DATA I DSC / 4*1 f 6*0 / 

DATA D6E0M / 10*0.0 / 

DATA JJ / 30*0 / 

END 

SubrOuYINe 

DUMMY SUBROUTINE TO SUPPRESS PRINTOUT OF TITLE PAGE. 

IF YOU DESIRE THIS PRINT f REMOVE THIS DUMMY SUBROUTINE. 
RETURN 
END 

"SUBROUTINE CONTES — 

DUMMY SUBROUTINE CALL BY 'LINKCAL ' ENTRIES IN INPUT DATA. 
RETURN 
END 

SUBROUTINE OUTVEC ( NNf ARRAY f TITLE ) 

THIS IS A SPECIAL OUTVEC FOR THE JUNCTURE CORNER ONLY. 
FOR THE MORE GENERAL CASEf REMOVE THIS SUBROUTINE. 

DIMENSION ARRAY < 1 ) f TITLE (8) 

CALL JNCVEC ( NNf ARRAY t TITLE ) 

RETURN 

END 

SUBROUTINE TBLINP 

* JNCINP ' INITIALIZES THE DATA FOR A JUNCTURE REGION. 

CALL JNCINP 

RETURN 

END 


APPENDIX A. 

SUBROUTINES FOR 3DPNS 
TEST CASE 


i«l 


0t ^Ous: quality 


jnpui data deck 


Vine FUSELAGE JUNCTURE (LOW. 

FENAME 
I NAME 0) 

NIZS 250, NSNUDI. > IV* MOSLEM 


ndp < io. li> » 


% 


NUAR •• it KMT PAS ■< 700, MKOtJt « 2. MUM7VC <> 10 »”h 0UI0* . l, 
NEO « V. MC0AV2 •< i, NEGAV3 - 1. NtiNKNN „ !i, 0 .. 4 } 

MPV8XT A NTABPT < 14 » NGMADD U. NEltJnl, 


\ JjJJJ 11 :J/ ‘ 6' LPPMCH - 2, JPR - i. NMBOUI • 50, 


\ MNROW 

> TMKQDS « 2* K0D6 


1V3» LCOE u 73, KNOW •« 75, WPVSX „ 2640. 

9r NRAND «■ 23, NODES '■ 1 70, NULL 14<>, 

it KODG « 1* KDUMP « 1* |<PN7 « i, 1WR1I L > . 
U-Nri^ 1 ^ 17 2 ’ IBOt U I,,)r W 4 * MR ‘ ■ i0 ’ NIT " 1 * NDBKP7 • G, 

ENAME02 

UINF - 100,0. F’INF « 2114,0. CHIEFS <« 3. Of. 4. HELP .. tol.O. 
V££*“ 2;» CO MPY « 1,. 1180 - 1,. EFMUET « 1, , DSTAKI ! l! 

U2STRS - 1.0* C4EDUU » .001, PCFAC7 « i,. WACl l, , 

tS HUL 1 ah 1 '« VL.DKL1 - 1.0. 00US0 l.OC-S* EPSHIN 1.011-5, 

TO .011. TD - , 50V. ISADEi « ,05. DELMLT « 1.00 *4, 

RARRAY ( 384 1 * 1,0. RARRAY ( 305 > .. 1.0, 

SEND 

tr.?,^ I ESTABLISH VECTOR ADDRESSES FROM DATA INPUT « 

lr } T INPUI DEPENDENT VARIABLE NUMBERS , 

1562378900. 10*0. 10*11 1 | 

LINK4 9 7 CALL POTENT TO ESTABLISH VECTORS EOR TYPE n r, r rr 

'•SIDD T ESTABLISH DIRECTION OF DIAGONALS. 

1 "2 T 

LINK2 14 T CALL D8CRT2 TO GENERATE DISCRETIZATION. 

E 9 9 T ' N0 * UK r - LEMtN ' 18 IN E1A »»• FOR EACH SUPER ILL ME NT. 

NEFS 9 9 T i N0 ' ° r KLKKEN, ° IW r 'XR. FOR EACH SUPER ELEMENT , 

Sm l*4 T 1 ' m UUF>ER KL EMC NTS. 4 UUADRILAT ERAL 

SELCN T SUPER ELEMENT NODE CONNECTIONS. 

9 10 3 1 12 13 19 2, 

1 3 14 IS IV 1& 17 10 7 

DEPVAR 28V 290 1248 7 VAR. TO BE DISTRIBUTED. 209=2, 2V0==Y, 1240*01 
0.00. 0.869545 .01 VIS .0 .025 .05 .025 .0 .OVVV ,0VVV 40 Ul ' 

• 2 1,33 .869565 . .01910 0.0, 1.15 1,35 1.15 1,35 

0.0 0.0 .01918 .05 0.0 0.025 .05 .075 0.0 .01918 
.0 .0009 ,01918 0 . OVVV O.OVVV .00 ,000V .09 ,0009 
0.0 0.0 .758 0.0 0,0 .5 .5 0.0 0,0 .720 .50 .703 0,0 

.783 0.0 ,58 .720 0.0 .71 I 

DONE T END OF INPUT FOR DSCRT 2 , 

RErtI ' . ^ * 3 “ 2A '* REARRANGE LOWER RIGHT CORNER ELEMENTS. 

REAI, i8i 171 "•a'.* 24 m M iw 1 SFT* m ' H Lln C0UME * kLtMEM,a * 
REA 110 100 90 S I 63 26 ' }n ’ fLtMEMT 10 interior corner. 

I ARRAY 14 325 7 SET HELEN » 325 

DESCRIPT 204 T THEE FOR OUTPUT SECTION. 

WIND / FU8ELASE JUNCTURE FLOW, 


DONE 


T END OF THEE 


APPENDIX B, 

DATA DECK FOR 3DPNS JUNCTURE 


I Ml ’U 


PA I A PI.CK 


r.i.i : ; - 

r.r i'M '• 


1 

t ' 1 « 


cow i tie i 

WIND / HISELACE 
PONE 

pescript 3 3 ? 

REFERENCE 
LLNCTH. *«*«..«• 

VELOCITY 

DENSITY 

riiKPLRAIURC. . . . 

i:M I HAI.PY 

IR07. SPEC. HEAT. 

VISCOSITY 

LOCAL PRESSURE. 
LOCAL SOLUTION. 

12PF1X 

NWCLOM H'S 

Xl/LREF 

DONE 

MPAFIA -1 

3*2. 2 2 162 


SCALAR IK l N 1 00 1 
JONCIOKIi I LOU. 


FOR 00 1 10 1 CLC I (ON. 


I JOPAR PARAML I HR MILES FOR Dili Fill. 


LiSH M 

INCUCH-1N 

H *K-b 

c *u -u 

• Hi M 1 M * M 

4 IN 

* M **•*««* « 4 « 


* f T/tJ 4 4 4 4 4 (4 

4 M 4 A 1 4 4 4 4 4 4 4 

*H/U. * 

i CM/U * , * , 4 * 

(UK/rU, 4 4 4 

4 M 4 ft 4 4 4 4 44*4 

1’ C/Mi& 1 * 4 4 4 4 4 


• RANK INI- 4 4 4 4 

4 N 4 6 4 4 * « 4 4 4 4 

• KELVIN 


• klTU/UBM* 4 • « 

4 4 4 4 4 4 4 

4 KJ/KL * 1 4 * 4 * 


.D fU/LBK •«. . 

4 N 4 ft 44444*4* 

4 KJ/KC |>M 1 4 

*N(ft( *4 *4* * 

.LBK/f T * 8 . . . 

« N 4 ft 4 4 4 4 4 4 * 4 

*MI S/Ml* 4 * . 4 

• Him 4 *44* 


« KOI 4 4 4 * 4 4 4 

KN*I/M:*4 *44*4 

* 1 Om 4 4 4 *.* 

4 MftCH NO* 4 4 « 

4 4 4 4 4 4 4 

4 li MUNCY (4*44 

4 INI 4 MftR* 4 

4 H21 4 4 4 4 4 4 *4 

* C22 4 * 

• Oltt 

4 1* 1 ( 4 4 4 4444 

♦ H.5 1*44444*4 

* 44*41**4 

• (.33 .. ...... 

* b ii « 1 4 | * M I 

• PXl/LKKP 4 * 4 

<EPU1U)M« 44 4 

.DXlli/LRIF , . 

4 Klf.ru fcE NO 


164 163* 3*2 164 163, 

3*2 165 2, 2 -175 3*2, 3*2 176 

2 2 169 168 167, 3*2 108 2, 5*2, 

X0NUH8 -1 

999, 5*200, 999, 200 4*43, 200 27 200 27 
200 10 200 10 10, 200 58 200 58 200, 

200 97 200 97 200, 200 30 200 30 200, 

200 38 200 38 38, 999, 39 4*36, 300 154 100 

398 4*11 186, 200 4*11 139, 11 12 14 85 47 


3*2 170 
3*2 177 
5*2, 5*2 


174, 
178, 
T 


;*7, 


135 

1 


122 , 


l&ORD T COUNTER-CLOCKWISE LIST Of BOUNDARY NODES. 

10*11 1, 17*110 20, 10*1-1 190, 18*1-10 171 T 

1 T FIXED NODES FOR DEP. VAR. Nil. 1 
DONE 

1 T FIX MALL NODES. 

T FIXED NODES FOR PEP. VAR. NO. 2 
DONE 

1 T FIX WALL NODES. 

T FIXED NODES FOR PEP. VAR. NO. 3 
DONE 

1 T FIX WALL NODES. 

T FIXED NODES FOR DEP. VAR. NO. 7 
DONE 

T FIX WALL ONLY FOR PP 
T FIXED NODES FOR DEP. VAR. NO. S 
DONE 


KBNO 
ADD 

19*110 
KBNO 2 

ADD 

19*110 
KBNO 3 

ADD 

19*110 
KBNO 7 

ADD 

19*110 
KBNO 
ADD 

19*110 
KBNO 
ADD 

9*110 10, 


1 

8 

10 

9 


T FIXED NODES FOR 
DONE " 
9*110 110 T 


DEP. VAR. NO. 9 


90 


APPENDIX B, 

DATA DECK FOR 3DPNS JUNCTURE 






INPUT DATA DECK 


C "miO T T 8FEClm lV SET8 0F 10 HOUW EACH, 

CN 190*II 1 \ SF ' ECIFY ALL M0DE H08 ‘ 1N THE »* SFJS. 

I ARRAY 47 19 T SET LG - 47 

LINKS 4 T CALL HI MEN TO N0N-DIKtN8I0NAL17E INPUT. 

■Sy 3 4l ' , C T Km! fT"' lt "°’ HS * H " ““ s - 

CN 19 T8 19 T T SET UP ™° VEC10RS 0F LEN °f« 19 FOR NUDPCP. 

CNTNDS 1 SPECIFY NODES IN ROTH VECTORS. 

10*1-10 100. 9*1-1 9, 10*110 100. 9*1-1 189 f 
1 ARRAY 47 2 T SET LG * 2 


LINK2 

23 T 

CALL CPGTUP TO READ 

IN UP DATA 



0.0 

0.00999 

0*02784 

0.04731 

0*07165 

0.10210 

0.14010 

0*19000 

0.46940 

0.01914 

0*46940 

0.42770 

0*38980 

0.36030 

0.33690 

0*31740 

0.38290 

0.03126 

0*38290 

0.35650 

0*32130 

0.29170 

0.26730 

0*24680 

0.29870 

0.04628 

0*29870 

0.28500 

0*25710 

0.22970 

0.20600 

0*18560 

0.21490 

0.06409 

0*21490 

0.21020 

0*19090 

0.16830 

0.14700 

0*12800 

0.13190 

0.08458 

0*13190 

0.13350 

0*12220 

0.10560 

0.08826 

0*07195 

0.05079 

0.13310 

0*05079 

0.05659 

0*05191 

0.04182 

0.02944 

0*01691 

-0.10000 

0.16090 

*0*10300 

-0.09411 

-0*08947 

0.08799 

-0.08917 

-0*09173 

-0.1788 

0.22240 

-0*1788 

-0.1664 

-0*1584 

-0.1522 

“0.1480 

0*14310 

-0.3151 

0.2908 

-0*3151 

-0.3012 

-0*2879 

-0.2739 

-0.2601 

-*0*2466 

-0.4296 

0.40250 

-0*4296 

-0.4155 

*0*3982 

•0.3783 

-0.3571 

0*3343 

-0.54460 

0.50 

-0*54460 

-0.53120 

-0*51030 

-0.48480 

-0.45640 

-0*42300 

-0.57430 

0.5975 

-0*57430 

-0.56140 

-0*33960 

-0.51270 

-0.48250 

-0*44880 

-0.54460 

0.71 

-0*54460 

-0.53120 

-0*51030 

-0.48480 

-0.45640 

-0*42500 

-0.4296 

T 

-0*4296 

-0.4155 

-0*3982 

-0.3783 

-0.357 

-0*3345 

LINKS 

30 T CALL 

NODFCP TO 

ESTABLISH 

CP TABLES 

♦ 



APPENDIX B. 

DATA DECK FOP 3DPNS JUNCTURE 


oimgsmai. is 

OF POOR QUALITY 


INPUT DAT A DECK 


T RESET CMTPT6 FOR SOLUTION. 

T RESET CHTNDS FOR SOLUTION, 

1 T 

47 19 T SET LC «• 19 

U T T CA ckll! 0 SSKVJ 0 8t ' ^ EUMEMT ,WU ' ^TACHED TO NODES. 
10 ?! C ? LL . T8 . L lrLI° 6 . l8,RI f. U . , . fc . U1 AL0WU ^LUNMS OF NODES. 


f CALL N0DPPR 8 RtOR'EL DLBUU. 

k 0° 1 T u?SS*\\ TUKM “»»• 

CALL WLFLXS 

1.0 T SET UEDGE «• 1.0 
CALL TRBTHK 
CALL WLFLXS 
CALL SETDIF 

T TITLES FOR OUT RUT DEPENDENT VARIABLES. 


U2/UKEF 
PP / PS1A0 
V'W' 


U3/UREF 

U' 

W' 


NU/NUREF 

U'V' 


TKE/TKEREF 

U'U' 


1248 2248 3248 1247 3248» 
6271 7271 8271 T 


SPECIFY PRINTOUT VARIABLES. 

8 1247 3248* 6248 7248 3271 4271 S271 , 


9 * 1 ' S 13*1 F V ULT1PLIERS AN ° 0FERA1ORS F0 * W«IHT VARIABLES. 
15 T CALL TRBTHK 

3 T CALL WLFLXS 
6 T CALL SETDIF 

T ORDER OF CALLS AT END OF STEP IN IMPLCT. 

T CON- 3 DUMMY CALL. 

95 0 . 9fc-3 T SET EPTEST 

4 T CALL RMLDST 

107 0 T SET ME1E2 - 0 

5 T FIXED NODES FOR DEP. VAR. NO. 5 

DONE 

a ’ T 1V ?tv?., iJLrrl FIX UALL AWl * FREE STREAM NODES. 

6 T FIXED NODES FOR DEP. VAR. NO. A 

DONE 

"t FmtsmLM ' mKs ‘ 

T EXIT PROGRAM 


APPENDIX B. 

DATA DECK FOR 3DPNS JUNCTURE 




1 1 £ unfcf fowls H y VkC.tnfrv Son.e typical entries an* tecorOed. 


p« r r*. r* „ *.»* " u~ ‘ 4 * - 

r„ r-, r *?,.'« ,v 

~ c- - fr « ~ 

c »■ ^ *c > r 

^ *r f !t n n"* ^ 

1 * 1*1 

c c o o c o a c q o n o 

-C (► O H f I K| f / *3 & *> : 


i 6 . . m « .* fi ^ rn ’ 

w: f« i> «-» * f*i f*3 o 

* © Ci V * * ♦ < 

" »■, f T T ♦ *3 V” 

I * t • . * I 


. q> ci * j 

HI," *T 0* < © (C O W f ♦ « • *1 * ^ ® J 


I f»» _ 

* cl^D 


o f- I* • . h *- *3 N *. K D (Mil O t! O O w r 
HflO 03 Ci 3 

w aj f <► *♦ <► ©* © o % c< © < 2 _ 

-*■ isss^ns? ,, ,3 
”i?.^iss?s?§ 2 ss?s?f?? 


o*vrf *»ftflcottOfinf;^f oon 
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